AD-A227  661 


GENERAL  INSTRUCTIONS  FOR  COMPLETING  SF  298 

The  Report  Documentation  Page  (RDP)  is  used  in  announcing  and  cataloging  reports.  It  is  important 
that  this  information  be  consistent  with  the  rest  of  the  report,  particularly  the  cover  and  title  page. 
Instructions  for  filling  in  each  block  of  the  form  follow.  It  is  important  to  stay  within  the  lines  to  meet 
optical  scanning  requirements. 


Block  1. 


Block  2.  Report  Date.  Full  publication  date 
including  day,  month,  and  year,  if  available  (e.g. 
1  Jan  88).  Must  cite  at  least  the  year. 

Block  3.  Type  of  Report  and  Dates  Covered, 
State  whether  report  is  interim,  final,  etc.  If 
applicable,  enter  inclusive  report  dates  (e.g.  10 
Jun  87  -  30  Jun  88). 

Block  4.  Title  and  Subtitle.  A  title  is  taken  from 
the  part  of  the  report  that  provides  the  most 
meaningful  and  complete  information.  When  a 
report  is  prepared  in  more  than  one  volume, 
repeat  the  primary  title,  add  volume  number, 
and  include  subtitle  for  the  specific  volume.  On 
classified  documents  enter  the  title 
classification  in  parentheses. 

Block  5.  Funding  Numbers.  To  include  contract 
and  grant  numbers;  may  include  program 
element  number(s),  project  number(s),  task 
number(s),  and  work  unit  number(s).  use  the 
following  labels: 


Contract 

Grant 

Program 

Element 


PR  -  Project 
TA  -  Task 
WU  -  Work  Unit 

Accession  No. 


Block  6.  AuthorfsL  Name(s)  of  person(s) 
responsible  for  writing  the  report,  performing 
the  research,  or  credited  with  the  content  of  the 
report.  If  editor  or  compiler,  this  should  follow 
the  name(s). 


Block  7. 


^forming  Qrqamzat 
l  Self-explanatory. 


Block  8.  Performing  Organization  Report 
Number.  Enter  the  unique  alphanumeric  report 
number(s)  assigned  by  the  organization 
performing  the  report. 


Block  9.  Soonsoring/Monitc 
Namesfsl  and  Addressees). 

Block  10.  Sponsorinq/Moni 
Report  Number.  (If  known) 


Self-explanatory. 


Block  11.  Supplementary.  Notes,  Enter 
information  not  included  elsewhere  such  as: 
Prepared  in  cooperation  with...;  Trans,  of ....  To 
be  published  in  ....  When  a  report  is  revised, 
include  a  statement  whether  the  new  report 
supersedes  or  supplements  the  older  report. 


Block  12a.  Distribution/Availablitv  Statement 
Denote  public  availability  or  limitation.  Cite 
any  availability  to  the  public.  Enter  additional 
limitations  or  special  markings  in  all  capitals 
(e.g.  NOFORN,  REL,  ITAR) 


DOD 


DOE 

NASA 

NTIS 


See  DoDD  5230.24,  "Distribution 
Statements  on  Technical 
Documents." 

See  authorities 

See  Handbook  NHB  2200.2. 

Leave  blank. 


Block  12b. 


DOD 

DOE 


NASA 

NTIS 


DOD  -  Leave  blank 

DOE  -  Enter  DOE  distribution  categoric' 
from  the  Standard  Distribution  for 
Unclassified  Scientific  and  Technical 
Reports 

NASA  -  Leave  blank 
NTIS  -  Leave  blank. 


Block  13.  Abstract.  Include  a  brief  (Maximum 
200  words)  factual  summary  of  the  most 
significant  information  contained  in  the  report. 

Block  14.  Subject  Terms.  Keywords  or  phrases 
identifying  major  subjects  in  the  report. 


Block  15.  Numbei 
number  of  pages. 

Block  16.  Price  C< 
code  (NTIS  only). 


Enter  the  total 


Enter  appropriate  price 


Blocks  17.-19.  Security  Classifications. 
Self-explanatory.  Enter  U.S.  Security 
Classification  in  accordance  with  U.S.  Security 
Regulations  (i.e.,  UNCLASSIFIED).  If  form 
contains  classified  information,  stamp 
classification  on  the  top  and  bottom  of  the  page. 

Block  20.  Limitation  of  Abstract.  This  block 
must  be  completed  to  assign  a  limitation  to  the 
abstract.  Enter  either  UL  (unlimited)  or  SAR 
(same  as  report).  An  entry  in  this  block  is 
necessary  if  the  abstract  is  to  be  limited.  If 
blank,  the  abstract  is  assumed  to  be  unlimited. 

Standard  Form  298  Back  (Rev.  2-89) 


A  FINITE  DIFFERENCE  NUMERICAL  ANALYSIS  OF  HEAT 
TRANSFER  IN  ATHEROMATOUS  PLAQUE  FOR 
PERCUTANEOUS  TRANSLUMINAL  MICROWAVE 
ANGIOPLASTY 


by 

Lori  Ann  Yuk  Lan  Young 


A  thesis  submitted  to  the  faculty  of 
The  University  of  Utah 

in  partial  fulfillaent  of  the  requirements  for  the  degree  of 


Master  of  Science 


Departaent  of  Mechanical 
The  University  of 


Engineering 

Utah 


90  10  9/A  08  8 


June  1990 


Copyright  ©  Lori  Ann  Yuk  Lan  Young  1990 
All  Rights  Reserved 


THE  UNIVERSITY  OF  UTAH  GRADUATE  SCHOOL 


SUPERVISORY  COMMITTEE  APPROVAL 


of  a  thesis  submitted  by 


Lori  Ann  Yuk  Lan  Young 


This  thesis  has  been  read  by  each  member  of  the  following  supervisory  committee  and  by  majority  vote 
has  been  found  to  be  satisfactory. 


Kuan  Chen 


THE  UNIVERSITY  OF  UTAH  GRADUATE  SCHOOL 


FINAL  READING  APPROVAL 


To  the  Graduate  Council  of  The  University  of  Utah: 


I  have  read  the  thesis  of  Lori  Ann  Yuk  Lan  Young _ in  its 

final  form  and  have  found  that  (1)  its  format,  citations,  and  bibliographic  style  are 
consistent  and  acceptable;  (2)  its  illustrative  materials  including  figures,  tables,  and 
charts  are  in  place;  and  (3)  the  final  manuscript  is  satisfactory  to  the  Supervisory 
Committee  and  is  ready  for  submission  to  the  Graduate  School. 


Dale 


Robert  F.  Boehm 

Member.  Supervisory  Committee 


Approved  for  the  Major  Depanmen; 


David  W.  Hoeppnef 

Chairman  Dean 


Approved  for  the  Graduate  Council 


"2 


B.  Gale  Dick 

Dean  of  The  Graduate  School 


ABSTRACT 


-4 

Various  percutaneous  transluminal  angioplasty  methods  have  been 
under  developaent  for  the  treatment  of  atherosclerosis.  The  microwave 
angioplasty  method  is  a  new  technique  that  has  a  treaendous  potential 
for  selectively  heating  atheroaatous  plaques  in  arteries  with  ainiaal 
thermal  damage.  The  purpose  of  this  study  is  to  evaluate  the  heat 
transfer  within  the  vascular  systea  while  being  exposed  to  a  high 
energy  electromagnetic  field  induced  by  a  microwave  applicator. 

A  heat  transfer  analysis  was  accomplished  using  the  explicit 
finite  difference  aethod.  The  angioplasty  system  was  modelled  as  a 
two-dimensional  cylindrical  coordinate  thermal  system  to  solve  for  the 
conductive  heat  transfer  in  the  atheroaatous  plaque,  vascular  tissue, 
and  muscle^  As  a  first  approximation,  the  convective  heat  transfer  was 
modelled  using  a  one-dimensional  cylindrical  coordinate  system  assuming 
steady,  laminar,  fully  developed  fluid  flow  and  Newtonian  fluid 
properties.  Three  models  with  different  microwave  applicator  designs 
were  analyzed  in  this  study.  All  three  models  were  analyzed  for  one 
second  of  heat  generation.  Further  studies  were  also  conducted  to 
analyze  the  effect  of  pulsing  the  microwave  heat  source. 

Results  of  this  study  confirm  the  theory  of  selective 
atheromatous  plaque  heating.  Two  of  the  three  models  produced 
temperatures  within  the  100  to  300 *C  range  to  induce  plaque  ablation. 

As  heat  dissipated  from  the  plaque  into  the  blood  vessel  and  muscle 
tissue,  some  thermal  damage  was  incurred.  Average  temperatures  in  the 


blood  vessels  (60  to  70*C)  were  sufficiently  high  to  initiate  tissue 
dehydration  and  shrinkage.  Hemolysis  due  to  a  significant  aaount  of 
heating  was  also  observed  in  the  blood.  In  one  of  the  models,  blood 
temperatures  as  high  as  400 *C  were  computed  in  the  thermal  anaylsis. 
Pulsing  of  the  microwave  heat  source  did  not  initially  reduce  blood 
vessel  wall  temperatures.  Since  the  blood  temperatures  were  often 
higher  than  that  of  plaque,  heat  was  always  conducted  in  the  outward 
radial  direction.  Consequently,  the  intima  layer  next  to  the  heated 
plaque  was  highly  susceptible  to  thermal  injury. 

The  future  of  microwave  angioplasty  is  promising.  Further 
development  of  the  applicator  design  may  eliminate  some  of  the 
hemolysis  and  tissue  damage  observed  in  this  preliminary  study.  The 
shape  of  the  microwave  applicator  had  a  significant  effect  on  the  heat 
generated  within  the  plaque  and  surrounding  tissue.  Development  of  a 
heat  exchanger  to  limit  adverse  blood  temperature  increases  during 
treatment  may  be  of  benefit. 
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CHAPTER  1 


INTRODUCTION 

Historical  and  General  Background 
The  use  of  heat  for  therapeutic  applications  has  been  exploited 
since  ancient  times.  For  exaaple,  heated  aetals  and  lenses  were  used 
to  control  the  size  of  surface  lesions  and  saall  tuaors,  and  treat 
other  diseases.  More  sophisticated  aethods  for  heating  tissues,  such 
as  lasers,  aicrowave  diathermy  and  ultrasound,  were  invented  during  the 
20th  century  and  continue  to  be  further  developed. 

Hypertheraia  treatment  is  an  application  of  aicrowave  diatheray 
in  which  a  aicrowave  and  radio-frequency  heat  generator  effectively 
heats  tumors,  both  malignant  and  nonmalignant,  to  remove  abnormal 
tissues.  Microwave  diatheray  involves  the  absorption  of  high-energy 
electromagnetic  waves  by  biological  tissues.  The  energy  carried  by 
these  waves  is  converted  into  thermal  energy  by  dielectric  and 
resistive  losses  (Jain,  1985).  By  elevating  the  temperatures  of  the 
tuaors  to  40-42 *C,  researchers  found  that  local  heating  of  the  tissues 
will  result  in  a  total  regression  of  melanomas  and  sarcomas  while 
increasing  the  survival  rate  of  patients.  The  success  of  this 
application  has  been  limited  by  the  superior  thermoregulation  system 
within  the  human  body,  thermal  properties  of  its  tissues,  and  thermal 
tolerance  of  the  tissues  to  withstand  the  higher  temperatures  without 


being  damaged. 
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Recent  studies  have  shown  that  the  majority  of  the  heat 
deposition  produced  by  the  electromagnetic  field  was  found  in  the  fat 
layer  of  the  body  (Christensen  and  Durney,  1981).  Although  this 
phenomenon  has  been  detrimental  in  the  application  of  hyperthermia  to 
shrink  tumors,  it  may  be  applied  to  treating  patients  with 
atherosclerosis.  If  the  atheromatous  plague  layer  absorbs  most  of  the 
heat  produced  by  the  microwave  applicator,  the  plague  cm  be 
selectively  ablated  without  causing  thermal  injury  to  the  blood  vessel 
and  surrounding  muscle  tissues. 

Atherosclerosis  is  a  condition  characterized  by  the  degeneration 
and  hardening  of  the  walls  of  arteries  and  valves  of  the  heart.  The 
disease  is  caused  by  an  accumulation  of  fatty  and  other  substances  in 
the  inner  lining  of  the  arteries  to  form  atheromatous  plagues  as  shown 
in  Figure  1,  These  plagues  intrude  upon  the  passageway  and  gradually 
obstruct  the  flow  of  blood  (Miller  and  Keane,  1987).  There  are  three 
basic  types  of  atheromatous  plagues:  (1)  calcified  plagues  that  are 
hard  and  bony,  (2)  fibrous  plagues  that  are  firm  and  pliable,  and  (3) 
fatty  plagues  that  are  loose  and  buttery  (Welch  et  al.,  1985).  This 
condition  is  a  major  cause  of  heart  disease  such  as  angina  pectoris, 
coronary  occlusion  or  coronary  thrombosis  (heart  attack),  and  cerebral 
vascular  accident  (stroke). 

Atherosclerosis  can  be  treated  by  in  vivo  percutaneous 
transluminal  angioplasty.  In  this  procedure,  a  vascular  catheter  with 
an  uninflated  balloon  is  inserted  into  the  femoral  artery  until  it 
reaches  the  point  of  blockage.  The  balloon  is  then  inflated  to 
mechanically  crush  the  plague  and  recanalize  the  blood  vessel.  The 
transluminal  angioplasty  method  for  treating  atherosclerosis  is 
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Figure  1.  Typical  cross  section  of  a  diseased  blood  vessel. 

preferred  over  the  coronary  bypass  in  terms  of  patient  recovery  and 
lowered  surgical  complications  as  well  as  reduced  medical  expenses 
associated  with  this  procedure. 

Modern  angioplasty  techniques  are  the  by-products  of  two 
important  technological  groundworks:  (1)  medical  imaging  based  on 
ionizing  energy  or  X-rays  and  (2)  access  into  the  vascular  lumen  using 
a  catheter.  The  X-ray  beam  was  first  discovered  by  Wilhelm  Konrad 
Roentgen  in  1895.  However,  it  was  not  until  1929  that  R.  Dos  Santos  et 
al.  demonstrated  the  use  of  X-rays  to  visualize  the  aorta  by  using  a 
contrast  material  injected  into  the  lumbar  aorta  (Jang,  1986).  Robb 
and  Steinberg  showed  the  heart  and  its  central  circulatory  system  could 
also  be  visualized  by  1938.  Several  major  advances  took  place  in  the 
1940s  with  Farinas  being  the  first  to  develop  an  aortogras  technique 
using  a  catheter  which  was  introduced  into  the  femoral  artery  and 
threaded  through  the  abdominal  aorta  (Jang,  1986).  By  1948,  a  radial 
catheter  technique  was  introduced  by  Radner.  Although  his  technique 
lacked  the  sophistication  and  quality  required  to  perform  transluminal 
angioplasty,  with  the  advent  of  TV  monitoring  systems  and  video 


recorders  to  provide  instant  replay  capability,  the  technology 
developed  rapidly  in  the  1950s  to  1970s.  In  1953,  a  set hod  for  placing 
a  catheter  into  the  vascular  systea  by  surgical  exposure  of  a  blood 
vessel  or  percutaneous  technique  was  developed  by  S.I.  Seldinger  (Jang, 
1986).  The  concept  for  applying  the  percutaneous  technique  to  perform 
transluainal  angioplasties  was  introduced  by  Dotter  in  1964.  The 
first  in  vivo  transluainal  coronary  angioplasty  was  performed  by 
Gruentzig  in  Zurich,  Switzerland  in  1977.  This  procedure  was 
introduced  to  the  United  States  in  1978  by  Myler  and  Stertzer.  Since 
1977  the  technique  for  performing  transluminal  angioplasty  has  advanced 
significantly  in  the  development  of  steerable  guide-wire  system  and 
improved  balloon  catheters  (Jang,  1986). 

Thus  far,  transluminal  coronary  angioplasty  has  not  shown  a 
significant  reduction  in  mortality  or  improvement  in  left  ventricular 
function.  In  a  six  year  study  conducted  by  Simpfendorfer  and 
associates,  approximately  4,000  patients  treated  by  transluminal 
coronary  angioplasty  between  the  years  1981  and  1987  showed  an  average 
primary  success  rate  of  90  percent  (Simpfendorfer  et  al.,  1988). 

Despite  the  relatively  high  initial  success  rate,  six  months  following 
treatment  about  one  third  of  the  patients  with  single  vessel  disease 
had  recurring  atherosclerotic  lesions  (restenosis) requiring  retreatment 
or  coronary  bypass  as  shown  in  Figure  2.  Moreover,  most  successful 
angioplasty  procedures  resulted  in  varying  degrees  of  vessel  wall 
injuries  that  occurred  at  the  site  of  arterial  wall  dilation.  These 
injuries  included  intimal  splitting,  subintimal  dissection,  medial 
tears,  and  submedial  dissection  as  shown  in  Figure  3  (Duber  et  al., 
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Figure  2.  Restenosis  in  a  blood  vessel  six  eonths  following  treatment 
with  translueinal  coronary  balloon  angioplasty. 
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Figure  3.  Typical  blood  vessel  cross  section  following  translueinal 
coronary  balloon  angioplasty  with  vessel  wall  injuries. 


Within  the  past  several  years,  euch  research  has  been  done  to 
coebine  translueinal  coronary  angioplasty  with  the  use  of  a  laser  to 
ablate  (or  vaporize)  the  atheroeatous  plaque.  Laser  angioplasty 
provides  two  ieproveeents  over  balloon  angioplasty  by  opening  blood 
vessels  difficult  or  nearly  impossible  to  clear,  and  reeoving  the 
plaque  by  vaporization,  thus  reducing  the  probability  of  subsequent 
restenosis.  In  a  study  by  Sanborn  and  Cueberland  (Goldsmith,  1987),  34 
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combined  laser-balloon  angioplasty  procedures  were  analysed  following  a 
six  month  period.  They  found  26  of  the  30  vessels  recanalized  by  laser 
angioplasty  remained  well-cleared  with  minimal  or  no  restenosis. 
(Goldsmith,  1987)  Although  the  ablative  properties  of  the  laser 
angioplasty  makes  this  procedure  an  excellent  inhibitor  of  restenosis, 
it  does  have  a  major  problem  of  causing  vessel  wall  perforation. 

Sumnarv  of  Current  Angioplasty  Research 

The  laser  angioplasty  is  a  relatively  new  technology  and  much 
iesearch  has  been  recently  accomplished  in  search  of  a  more  reliable 
clinical  application  of  this  method.  For  example,  much  experimental 
work  has  been  done  to  characterize  various  types  of  laser  sources  such 
as  the  C02>  Nd-YAG,  and  argon  lasers.  In  1982,  Abela  et  al.  evaluated 
the  performance  of  these  three  laser  types  in  rabbits.  In  general, 
they  found  that  approximately  20  to  30  joules  of  energy  were  needed  to 
successfully  ablate  atherosclerotic  tissues.  They  were  also  able  to 
characterize  three  zones  of  thermal  damage:  (1)  complete  tissue 
vaporization,  (2)  thermal  burning,  and  (3)  diffuse  tissue  disruption. 
(Abela,  1982).  In  a  1987  thermographic  study  of  these  same  lasers, 
Nnitentag  et  al.  found  that  the  COj  laser  produced  the  best  heat 
absorption  while  the  Nd-YAG  laser  had  the  best  tissue  penetration 
characteristics  (Mnitentag  et  al.,  1987). 

To  minimize  the  complications  of  vascular  wall  perforation  and 
thermal  damage,  others  have  experimented  with  variations  to  the 
original  laser  angioplasty  technique.  For  example,  much  research  and 
experimentation  has  been  done  with  laser  heated  metal  tips  to  better 
control  the  laser  heating  within  the  blood  vessels.  In  an  attempt  to 
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understand  ablation  phenomena  Welch  et  al.  (1987)  studied  the  thermal 
characteristics  of  laser  probe  ablation.  Using  a  thermographic  and 
histologic  analysis  they  found  ablation  in  atheromatous  plaques  started 
at  100  *C  and  may  exceed  temperatures  up  to  180 'C  with  heating  of  the 
adventitial  tissue  to  temperatures  in  the  70’C  range.  In  addition, 
ablation  was  more  extensive  in  atheromatous  tissue  rather  than 
nonatherosclerotic  aortic  tissue.  Another  experimental  study  done  by 
Morcos,  Berns  and  Henry  revealed  that  the  laser  heated  tip  method 
increased  the  incidence  of  plaque  coagulation  and  vasospasm,  a  sudden 
but  transitory  constriction  of  the  blood  vessel  persisting  for  several 
hours  after  treatment  (Morcos  et  al.,  1988).  Another  variation  of  the 
laser  angioplasty  consists  of  combining  laser  heating  with  the  coronary 
balloon  method  to  mechanically  dilate  the  blood  vessel  while  thermally 
softening  the  plaque.  Preliminary  experimental  studies  conducted  by 
White  et  al.  demonstrated  that  this  method  produced  results  similar  to 
the  heated-tip  and  free  argon  laser  methods  and  had  no  significant 
advantage  over  other  laser  angioplasty  techniques  due  to  the  lack  of 
selective  heating  (White  et  al.,  1988).  In  a  1990  paper  published  by 
IEEE,  Rosen  proposed  a  similar  technique  for  combining  vascular 
dilation  and  thermal  ablation  using  a  microwave  heat  source  instead  of 
a  laser,  and  conducted  experiments  using  a  rabbit  model.  In  his 
application  of  this  technique,  the  microwave  heat  source  is  used  to 
soften  the  plaque  to  prevent  vascular  restenosis  rather  than 
selectively  ablate  the  atheromatous  plaque. 

While  many  attempts  have  been  made  to  develop  an  effective,  long 
term  treatment  for  atherosclerosis,  none  of  these  procedures  will 
selectively  ablate  atheromas.  A  better  understanding  of  the  thermal 
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and  physical  behavior  of  atheromatous  plaques  is  needed  to  devise  a 
■ethod  that  discriminates  biological  aaterial  types  and  removes 
abnormal  tissues.  Thus  far,  most  recent  research  efforts  have  been 
experimental  studies  characterizing  different  types  of  laser  sources 
for  fiber  optic  delivery  and  their  capability  for  plaque  ablation, 
while  studying  the  arterial  response  to  find  ways  of  controlling  the 
laser  to  minimize  the  damage  to  the  normal  tissues.  The  electrical  and 
thermal  properties  of  atheromatous  plaques  are  sketchy  at  best. 

Research  efforts  to  evaluate  the  thermal  properties  of  atheromatous 
plaques  have  been  limited  to  studies  conducted  by  Welch  et  al.  at  the 
University  of  Texas  in  1985  where  the  average  thermal  conductivity  and 
diffusivity  coefficients  of  various  types  of  plaques  were 
experimentally  determined  (Welch  et  al.,  1985). 

Many  numerical  studies  have  been  accomplished  to  study  the 
bioheat  equation  and  its  application  to  predict  the  heat  transfer 
effects  due  to  blood  perfusion  in  various  parts  of  the  body  and  in  the 
study  of  hyperthermia  to  treat  malignant  tumors.  However,  numerical 
modelling  of  angioplasty  procedures  to  predict  the  physical  behavior  in 
the  atheromatous  plaques  and  surrounding  vascular  tissues  has  been 
limited.  In  1987,  Nickolay  Furzikov  developed  a  quantitative  estimate 
for  predicting  tissue  destruction  based  on  laser  input  power,  time 
exposure,  radiation  spot  size,  tissue  thermal  conductivity,  and 
attenuation  coefficient.  His  thermal  model  agreed  with  experimental 
results  although  it  was  a  crude  approximation  due  to  the  lack  of 
reliable  data  for  the  physical  parameters  used  in  his  analysis 
(Furzikov,  1987).  Decker-Dunn,  at  the  University  of  Utah  in  1989, 
developed  a  mathematical  model  to  predict  normal  tissue  damage  during 
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laser  angioplasty  using  a  finite  difference  method.  The  Kubelka-Munk 
two-flux  theory  and  modified  Beer’s  Law  were  used  to  aodel  the  laser 
light  incident  froa  the  laser  tip  (Decker-Dunn,  1989).  To  date,  no 
nuaerical  theraal  analysis  has  been  done  for  the  aicrowave  angioplasty 
technique. 

This  project  investigates  the  use  of  a  aicrowave  heat  source  (as 
opposed  to  a  laser)  to  ablate  the  atheroaatous  plaque  in  blood  vessels. 
A  heat  transfer  analysis  is  accoaplished  using  a  finite  difference 
nuaerical  analysis  in  a  two-diaensional  cylindrical  coordinate  systea 
to  predict  the  teaperature  distribution  within  the  blood  vessel.  The 
result  of  this  theraal  analysis  is  a  first  approxiaation  characterizing 
the  theraal  behavior  of  a  aicrowave  percutaneous  transluainal 
angioplasty  systea.  An  analysis  is  coapleted  to  evaluate  atheroaatous 
plaque  ablation,  therapeutic  benefits,  possible  theraal  injuries  caused 
by  this  technique,  and  the  feasibility  of  applying  this  technique 
clinically  in  the  future. 

Further  developaent  of  the  finite  difference  heat  transfer 
analysis  is  needed  to  increase  its  accuracy  for  clinical  application. 
Refineaents  include  a  three-dimensional  solution,  detailed  analyses  of 
complex  theraal  and  fluid  phenomena  such  as  tissue  perfusion,  metabolic 
heat  generation,  and  non-Newtonian  blood  flow  effects  in  the  system, 
and  comparative  analysis  with  in  vitro  experimental  results.  With 
these  improvements,  cardiologists  may  use  this  thermal  analysis  to 
predict  teaperature  distributions  and  prescribe  the  therapeutic 
aicrowave  doses  required  to  achieve  adequate  and  safe  atheroaatous 
plaque  ablation  prior  to  surgical  procedures. 


CHAPTER  2 


HEAT  TRANSFER  THEORY 


Description  of  the  Angioplasty  Model 
The  physical  system  evaluated  consists  of  four  layers  of 
dissimilar  tissues,  which  include  the  blood,  atheromatous  plaque,  blood 
vessel,  and  muscle  representing  the  semi infinite  tissue  surrounding  the 
blood  vessel.  A  two-dimensional  cylindrical  coordinate  system  is 
applied  to  the  numerical  model.  As  an  initial  approximation,  the  model 
is  axially  symmetric  with  the  microwave  applicator  located  at  the 
center  of  the  blood  vessel.  The  microwave  applicator  is  used  with  a 
balloon  to  restrict  blood  flow  and  prevent  the  blood  from  heating  up 
too  rapidly  relative  to  the  atheromatous  plaque.  A  schematic  of  the 
model  analyzed  is  shown  in  Figure  4. 
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Figure  4.  Schematic  of  the  microwave  angioplasty  model 


In  1948,  Pennes  derived  the  bioheat  equation  to  analyze  the 
tissue  and  blood  temperatures  in  a  resting  human  forearm.  The  bioheat 
equation  is  widely  used  in  the  analysis  of  many  heat  transfer  problems 
with  biomedical  applications.  The  general  bioheat  equation  is  (Jain, 
1985): 

pcf*  =  7(k<7T)  +  qb  +  q.  (2-1) 

where:  T  =  tissue  temperature 

c  =  tissue  specific  heat  capacity 

p  =  tissue  density 

k  =  tissue  thermal  conductivity 

qb  =  rate  of  heat  exchange  with  the  blood 

q_  =  rate  of  metabolic  heat  generation 

t  =  time 

The  bioheat  equation  is  an  energy  balance  equation  about  a 
control  volume,  within  a  biological  system.  The  term  on  the  left  side 
of  the  equation  represents  the  transient  amount  of  energy  stored  by  the 
control  volume  which  produces  a  temperature  distribution  change  over  a 
period  of  time,  t.  The  first  term  on  the  right  side  of  the  equation 
accounts  for  the  heat  being  conducted  through  the  various  layers  of 
tissues  with  differing  thermal  properties  using  Fourier’s  Law.  For 
this  particular  model,  conductive  heat  transfer  is  calculated  for  the 
muscle,  blood  vessel,  and  atheromatous  plaque. 

The  second  term  of  the  bioheat  equation  accounts  for  the  heat 
transfer  due  to  the  blood  flow  (also  referred  to  as  perfusion)  within 
the  body’s  circulatory  system.  In  hyperthermia  applications,  the 
convective  heat  transfer  term  is  normally  used  to  account  for  the  net 
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thermal  energy  of  the  arterial  blood  flow  and  the  venous  blood  flow. 

For  the  angioplasty  aodel,  there  are  two  types  of  heat  transfer  caused 
by  the  blood  flow  to  the  region  being  evaluated:  (1)  blood  flow  within 
the  blood  vessel  and  (2)  tissue  perfusion.  Since  heat  transport  for  a 
single  blood  vessel  is  of  primary  interest,  this  analysis  is  liaited  to 
evaluating  the  heat  transfer  caused  by  the  arterial  blood  flow.  In 
this  case,  the  convective  heat  transfer  boundary  conditions  of  the 
blood  vessel  wall/atheroaatous  plaque  and  nicrowave  applicator  are 
luaped  into  the  governing  bioheat  equation  for  the  blood  flow.  The 
convective  heat  transfer  ter*  q^  can  be  expressed  in  teras  of  the 
convective  heat  transfer  coefficients  for  the  vascular  wall  and  plaque, 
hp,  and  aicrowave  applicator,  hj,  the  differential  peripheral  surface 
area  of  the  luaen,  dA,  and  the  blood  teaperature,  Tb,  surface  tissue 
teaperatare,  Tp,  and  aicrowave  applicator  surface  teaperature,  Tt ,  as 
shown  in  equation  (2-2). 

q„dV  =  hpdAp(Tb-Tp)+h.dA,(Tb-T.)  (2-2) 

Finally,  a  third  tera  is  used  to  represent  the  heat  that  is 
generated  due  to  natural  aetabolic  processes  in  the  body  and  external 
heat  generation  sources.  The  aetabolic  heat  generation  refers  to  the 
heat  produced  in  the  body  as  a  result  of  oxygen  consuaption.  This  fora 
of  heat  generation  is  calculated  by  aultiplying  the  oxygen  consuaption 
by  the  caloric  value  of  the  oxygen  (Jain,  1985).  In  this  theraal 
analysis,  calculation  of  the  teaperature  distribution  in  the 
atheroaatou8  plaque  (which  does  not  produce  any  aetabolic  heat 
generation)  is  of  priaary  interest.  Because  the  aaount  of  aetabolic 
heat  generated  in  the  theraal  aodel  is  relatively  saall  in  coaparison 
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to  the  heat  absorbed  by  the  tissues  due  the  electromagnetic  field 
induced  in  the  system,  it  be  treated  as  a  negligible  heat  source  and  is 
therefore  omitted  from  this  analysis. 

The  major  contributor  of  heat  generation  in  the  system  is 
produced  by  a  microwave  applicator.  The  surrounding  blood  and  tissues 
absorb  high-energy  electromagnetic  waves  produced  by  the  applicator, 
and  this  energy  is  converted  into  thermal  energy  by  the  dielectric  and 
resistive  losses.  The  amount  of  microwave  energy  generating  thermal 
energy  in  the  tissues  is  a  function  of  the  electrical  conductivity  of 
the  tissue,  o  and  the  electric  field,  E  (Jain,  1985). 

Q  =  JoE2  (2-3) 

The  electric  field  in  equation  (2-3)  is  computed  using  the 
Maxwell  equations  which  form  the  framework  of  classical  electromagnetic 
theory.  For  this  particular  application  the  finite  difference  time 
domain  (FDTD)  technique  was  used  to  solve  these  equations.  Heat 
generated  by  the  microwave  applicator  was  solved  directly  using 
equations  (2-4)  and  (2-5). 

VxE  =  -»g  (2-4) 

VxH  =  J  +  eg  (2-5) 

where:  E  =  electric  field 
H  =  magnetic  field 
U  =  material  permeability 
£  =  complex  permittivity 


t  =  time 


14 


In  general,  aicrowave  di at he ray  studies  applied  in  hyperthermia 
treataents  have  shown  that  tissues  containing  high  water  concentrations 
such  as  auscle,  skin,  and  internal  organ  tissues  have  high  absorption 
and  low  depth  of  penetration  characteristics.  On  the  other  hand, 
tissues  with  low  water  content  such  as  fat  and  bone  tissue  tend  to  have 
low  absorption  and  high  penetration  properties  (Jain,  1985). 

As  a  first  approxiaation  for  applying  the  bioheat  equation  to  the 
transluainal  aicrowave  angioplasty  aodel,  aodif ications  to  the  basic 
equations  were  required.  First  of  all,  the  heat  transfer  aodel  was 
reduced  to  a  two-diaensional  axially-syaaetric  cylindrical  coordinate 
system  (r,x)  reducing  the  bioheat  equation  to  the  following  fora. 


where:  q^  =  heat  generated  by  the  aicrowave 
applicator 

ot  =  tissue  theraal  diffusivity 
k  =  tissue  thermal  conductivity 
T  =  tissue  temperature 
t  =  time 

qb  =  heat  transfer  of  the  blood 


The  9  term  is  eliminated  from  the  governing  equation  since 
symmetry  with  respect  to  the  axial  direction  is  assumed,  and  the 
temperature  distribution  does  not  to  vary  significantly  in  this 
direction  of  the  theraal  aodel.  In  actual  situations  where  the 
geoaetry  of  the  atheromatous  plaque  aay  be  highly  irregular,  this 
assumption  aay  not  lend  itself  to  an  accurate  approximation.  However, 
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as  an  initial  analysis  to  predict  the  heat  distribution  of  the 
transluainal  Microwave  angioplasty  systea  and  evaluate  the  feasibility 
of  its  application,  a  two-diaensional  systea  will  be  analyzed. 

Finite  Difference  Solution  Of  the  Bioheat  Equation 

Because  of  the  highly  irregular  geoaetry  of  the  atheroaatous 
plaques,  an  analytical  solution  for  the  angioplasty  aodel  is  not 
available.  A  finite  difference  Method  is  applied  to  solve  nuaerically 
the  bioheat  equation.  There  are  three  Methods  for  evaluating  the  tiae 
derivative:  (1)  the  pure  iaplicit  Method,  (2)  the  Crank-Nicolson 
implicit  aethod,  and  (3)  Euler  explicit  Method.  The  iaplicit  Methods 
are  generally  More  advantageous  because  they  do  not  produce  the 
numerical  oscillations  that  may  occur  if  too  large  a  tiae  step  is  used. 
However,  the  explicit  Method  was  used  in  view  of  its  siaplicity  and 
fewer  number  of  calculations.  Unlike  the  iaplicit  Methods,  a  gauss i an 
simultaneous  equation  solver  subroutine  is  not  required  because  the  new 
temperatures  are  computed  directly  from  the  physical  foraulation 
equations  derived. 

To  apply  the  explicit  finite  difference  Method,  a  two  dimensional 
grid  was  established  with  semi infinite  and  convective  boundary 
conditions  as  shown  in  Figure  5. 

The  finite  difference  program  developed  for  this  thesis  project 
computes  the  transient  heat  and  mass  transport  only.  Another  finite 
difference  program  to  compute  the  energy  deposition  from  the  microwave 
applicator  developed  by  Paul  McArthur  and  Wenquan  Sui  of  the  University 
of  Utah,  Hyperthermia  Research  Group,  Department  of  Electrical 


16 


Senl- 

Inflnlte 

Boundary 


Hlcro«Mue 

Applicator 


Figure  5.  Axisymmetric  percutaneous  transluminal  microwave  angioplasty 
thermal  aodel  with  semiinfinite  boundary  conditions. 


Engineering  was  aerged  with  this  prograa  to  predict  the  overall 
temperature  distribution  of  the  microwave  angioplasty  aodel. 

The  thermal  models  analyzed  were  generated  by  a  aodel  generator 
prograa  developed  by  Paul  McArthur  of  the  Department  of  Electrical 
Engineering.  Modifications  to  the  coordinate  system  were  made  to 
accommodate  the  heat  deposition  data  passed  to  this  prograa  in  order  to 
compute  the  temperatures  at  the  boundary  interfaces  to  solve  accurately 
for  temperatures  where  both  conduction  and  convective  heat  transfer  was 
occurring  at  the  same  node.  This  was  accomplished  by  shifting  the  two 
dimensional  cylindrical  grid  system  by  DR/2  and  DX/2  in  the  x  and  r 
directions,  respectively.  The  new  heat  deposition  values  at  these  new 
points  were  calculated  using  the  cubic  spline  interpolation  method. 

The  numerical  solution  consists  of  twelve  different  types  of 
finite  different  nodes  to  simulate  the  conduction  heat  transfer  process 
within  the  muscle,  blood  vessel,  and  atheromatous  plaque.  Each  node 
type  was  derived  using  the  physical  formulation  method  to  account  for 
the  various  internal,  boundary,  and  multiproperty  nodes  modelled  in  the 
analysis.  Calculations  for  the  coefficients  of  the  temperature  at  the 
node  being  evaluated  (T-  =  )  were  also  done  using  the  Biot  and  Fourier 
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Numbers.  This  calculation  was  extremely  valuable  in  preventing  the 
occurence  of  unstable  oscillations  in  any  part  of  the  numerical 
solution. 

The  energy  balance  and  finite  difference  equations  for  each  of 

the  twelve  nodes  used  in  the  finite  difference  analysis  are  shown  in 

Figures  6  through  17  on  the  following  pages  The  first  set  of  physical 
formulations  (Figures  6  to  11)  represent  the  finite  difference  nodes 
occurring  at  the  insulated  and  interior  layers  of  tissue  in  the  thermal 
model.  The  nodes  displayed  in  Figures  12  to  17  nodes  are  used  at  the 
convective  boundary  layers  of  the  atheromatous  plaque  or  blood  vessel. 
At  the  boundary  interfaces,  semiinfinite  conditions  were  assumed  at 
both  ends  in  the  x  direction  and  for  relatively  large  radius  values 

where  the  muscle  tissues  are  not  affected  by  heat  conduction  and  the 

high-energy  electromagnetic  waves.  The  nodes  at  these  boundaries  are 
approximated  as  insulated  surfaces. 

In  figures  showing  more  than  one  configuration  the  energy  balance 
equation  is  given  for  the  figure  on  the  left  side.  The  other  equation 
may  be  easily  derived  by  changing  the  subscripts  for  the  symmetrical 
case.  For  the  finite  difference  cells  containing  one  material  type,  BO 
refers  to  the  inner  radius  closed  to  the  microwave  applicator  and  R1 
refers  to  the  outer  radius.  In  cells  with  two  material  types,  RO  is 
again  the  inner  radius,  R1  represents  the  radius  at  the  interface  of 
the  two  dissimilar  materials,  and  R2  refers  to  the  outer  radius. 
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6.  Interior  node  of  a  single  aediua  (Subroutine  INTRIOR). 
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Figure  7.  Interior  node  with  two  Media  (Subroutine  R2MATL) 


19 


or 


Ar 


TT-'.-n 

m 

i±Ul 

Zi 

m 


m 


V,. 


A  T  ’  -T,  ,  T  _t 

nrfvl.  2  ,  ZiAx,  '  •  i  ’ J 

Pc[fl(r,  -r0  )  2  ] 

k^(r,z-r02) 

i  rr  m  m 

ILli  .  1  1  •  i  ,  „  T 2  „  2V&X 

+  Ar  ♦  0,^(1*!  -r0  )  2 


_ Ar _ 

k^TjAx 


k*r0Ax 


(2-9) 


Figure  8.  Seui infinite  boundary  node  of  a  single  aediua  at  either  end 
of  the  aodel  in  the  x  direction  (Subroutine  SIDE1). 
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Figure  9.  Sesi infinite  boundary  node  with  two  Media  at  either  end  of 
the  Model  in  the  x  direction  (Subroutine  SIDE2). 
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Figure  10.  Sesi infinite  boundary  interface  node  in  the  auscle  tissue 
(Subroutine  INSBND). 
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Figure  11.  Sesi infinite  boundary  corner  node  in  the  uuscle  tissue 
(Subroutine  INSCNR). 
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Figure  12.  Convection  at  the  blood  vessel  corners  (Subroutine 
CONVCNR ) . 
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Figure  13.  Convection  at  the  blood  vessel/blood  or  plaque/blood 
interfaces  (Subroutine  CONVBND). 
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Figure  14.  Subroutine  CONV3A  or  CONV3B. 
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Figure  15.  Subroutine  CONV4A  or  C0NV4B. 
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Figure  16.  Subroutine  CONV5A  or  CONV5B 
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Figure  17.  Subroutine  C0NV6A  or  CONV6B 
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Heat  transfer  to  or  froa  the  blood  flow  was  calculated  in  a 
separate  subroutine.  As  an  initial  approxiaation,  the  convective  heat 
and  aass  transfer  of  the  blood  was  solved  nuaerically  using  a  one 
diaensional  finite  difference  aodel.  The  convective  heat  transfer 
coefficient  was  estiaated  using  the  two  analytical  solutions:  (1)  the 
solution  for  a  concentric  annular  pipe  which  gives  the  Nusselt  Nuabers 
as  a  function  of  the  ratio  of  the  aicrowave  applicator  radius  and  blood 
vessel  or  plaque  radius  and  (2)  the  solution  for  circular  pipe  flow 
(see  equation  2-20).  In  order  to  use  these  solutions,  the  teaperature 
and  velocity  profiles  of  the  blood  flow  were  assuaed  to  be  steady  and 
fully  developed  laainar  flows.  A  cubic  spline  interpolation  aethod  was 
used  to  calculate  the  average  Nusselt  nuabers  between  the  values  given 
in  Table  1  for  the  concentric  annular  pipe  solution.  The  convective 
heat  transfer  coefficient  is  then  calculated  knowing  the  hydraulic 
diaaeter  and  theraal  conductivity  by  the  following  equations. 


Nu  =  nr  =  3,66 

(2-20) 

Dh  =  2(r0-r, ) 

(2-21) 

Jjef 

ii 

J3 

(2-22) 

A  siailar  derivation  of  the  finite  difference  energy  balance 
equation  for  each  one-diaensional  eleaent  was  conducted  to  siaulate  the 
heat  transfer  contribution  due  to  the  blood  flow  in  the  systea.  The 
finite  difference  equation  accounts  for  the  advection  of  the  blood  flow 
through  each  control  voluae,  the  convective  heat  transfer  froa  the 
aicrowave  applicator  and  vascular /plaque  boundary  conditions,  and 
transient  heat  storage  of  the  blood.  Figure  18  shows  a  typical  one 
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Table  1.  Circular-Tube- Annulus  Convective  Heat  Transfer  Solution 


r* 

Nu, 

Nu, 

0.00 

« 

4.364 

0.05 

17.81 

4.792 

0.10 

11.91 

4.834 

0.20 

8.499 

4.883 

0.40 

6.583 

4.979 

0.60 

5.912 

5.099 

0.80 

5.58 

5.24 

1.00 

5.385 

5.385 

Note:  Analytical  solution  for  the  concentric  circular-tube  annulus 
with  laminar, fully  developed  velocity  and  tenperature  profiles  as  a 
function  of  r*=rB/r  where  the  subscript  a  refers  to  the  microwave 
applicator  and  p  refers  to  the  plague  or  vascular  wall.  Nu  and  Nu, 
represent  the  respective  Nusselt  numbers.  (Adapted  from  Kays  & 
Crawford,  1980) 
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Figure  18.  One-dimensional  finite  difference  node  for  modelling 
convective  heat  transfer  due  to  vascular  blood  flow. 
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dimensional  finite  difference  element  with  its  accompanying  governing 
energy  equation  derived  using  the  physical  formulation  method. 

The  maximum  possible  velocity  of  the  blood  was  computed  assuming 
an  average  flow  rate  of  250  ml/min  and  dividing  by  the  difference  of 
the  cross-sectional  area  of  the  lumen  and  the  area  of  the  microwave 
applicator.  Blood  mass  flow  rates  vary  depending  on  the  rate  the  heart 
is  pumping  and  the  part  of  the  body  where  the  blood  vessel  is  located. 
Since  this  program  may  be  used  for  applying  the  microwave  angioplasty 
technique  for  other  parts  of  the  body  with  lower  mass  flow  rates,  m, 
the  computer  code  has  been  written  to  accommodate  other  thermal  models 
with  slower  blood  flow  rates. 

As  a  first  approximation  the  thermal  model  of  the  blood  flow  does 
not  account  for  the  heat  transfer  caused  by  tissue  perfusion,  the  rate 
of  heat  exchange  occurring  in  the  microvasculature  of  the  surrounding 
muscle  tissues  (Shitzer,  1985).  Microcirculation  provides  the  means  by 
which  heat  generated  by  the  microwave  applicator  is  dissipated  within 
the  muscle  layers.  For  the  transluminal  angioplasty  model,  heat 
transfer  due  the  microvascular  blood  flow  is  negligible  relative  to  the 
heat  being  generated  in  the  system  and  conducted  to  the  blood  within 
the  vessel  during  treatment.  The  mathematical  models  using  the  Beer’s 
Law  and  Kubelka-Munk  two- flux  theory  developed  by  Decker-Dunn  for  the 
laser  angioplasty  model  also  confirms  this  assumption.  Her  studies 
showed  that  there  were  no  significant  theoretical  changes  by  modelling 
the  tissue  perfusion  contribution  in  her  thermal  analysis  (Decker-Dunn 
et  al. ,  1989) . 

In  modelling  the  blood  flow  it  was  assumed  the  flow  throughout 
the  system  was  laminar.  The  Reynolds  number  was  calculated  using  an 
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average  artery  diaaeter  of  0.4  ca,  total  blood  flow  rate  in  the  huaan 
body  of  250  al/ain,  blood  density  of  1.057  ga /cmP,  and  a  dynamic 
viscosity  of  0.034  ga/ca  sec.  This  yielded  a  value  for  the  Reynolds 
nuaber  of  412,  which  is  considerably  less  than  the  critical  Reynolds 
nuaber  of  2000  for  laainar  flow.  With  a  fifty  percent  occluded  artery 
(artery  diaaeter  of  0.2828  ca)  and  a  aicrowave  applicator  diaaeter  of 
0.2  ca,  using  the  hydraulic  diaaeter  for  a  concentric  annular  tube,  the 
calculated  Reynolds  nuaber  is  1991.  This  is  the  aaxiaua  possible 
Reynolds  nuaber  since  blood  flow  in  the  arteries  is  usually  less  than 
250  al/ain  and  varies  considerably  depending  upon  the  individual  and 
the  location  of  the  blood  vessel.  Although  blood  is  a  non-Newtonian 
fluid  since  the  viscosity  of  blood  varies  with  the  percentage  of  red 
blood  cells  present  (the  heaatocrit),  blood  viscosity  does  not  change 
significantly  over  nornal  ranges  of  flow  in  larger  blood  vessels 
(Selkurt,  1982).  Therefore,  for  the  purposes  of  this  study  a  constant 
viscosity  was  assuaed  and  the  blood  flow  was  treated  as  a  Newtonian 
fluid. 

The  explicit  finite  difference  equations  with  ablation  are 
siBilarly  derived.  See  Appendix  D  for  a  detailed  analysis  of  the 
development  of  these  equations. 


CHAPTER  3 


NUMERICAL  ANALYSIS  AND  RESULTS 

Investigative  Procedures  and  Graphical  Results 
The  Finite  Difference  Angioplasty  Thermal  Analysis  (FDATA) 
program  was  applied  to  three  different  microwave  angioplasty  models 
generated  by  the  University  of  Utah  Hyperthermia  Research  Group.  Heat 
generation  data  were  computed  using  a  Finite  Difference  Time  Doaain 
(FDTD)  Method  to  solve  for  the  electric  and  Magnetic  fields.  These 
data  were  read  into  the  FDATA  program  to  compute  the  temperature 
distribution  for  one  second  heat  generation  by  the  microwave 
applicator.  The  heat  generation  data  calculated  by  the  Hyperthermia 
Group  assume  steady  state  conditions.  Therefore,  the  heat  absorbed  by 
various  tissues  was  calculated  independently  of  the  FDATA  program.  The 
heat  generated  at  each  node  is  treated  as  a  step  function  where  heat  is 
suddenly  turned  on  and  maintained  at  a  constant  for  one  second. 

The  first  two  models  analyzed  are  pictured  in  Figures  19  and  20. 

A  200  by  30  cylindrical  grid  system  was  constructed  with  a  cell  size  of 
0.14286  cm  in  both  the  x  and  r  directions.  A  corresponding  time  step 
of  2.5  x  10'*  seconds  was  used  without  causing  any  unstable  oscillations 
in  the  finite  difference  solution.  Three-dimensional  graphs  of  the 
initial  specific  absorption  rate  or  SAR  (see  Appendix  A  for  a  more 
detailed  discussion  of  selective  microwave  heating)  distribution  for 
Models  1  and  2  are  shown  in  Figures  21  and  22,  respectively.  The 
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Figure  19.  Diagram  of  Model  1 


Figure  20.  Diagram  of  Model  2 

resulting  temperature  distributions  after  one  second  of  microwave 
heating  is  exhibited  in  Figures  23  through  26.  The  results  are 
displayed  as  both  three-dimensional  surface  and  two-dimensional  contour 
graphs. 

A  third  model  was  generated  to  simulate  better  actual  dimensions 
within  the  vascular  system.  The  model  was  generated  using  typical 
dimensions  for  am  average  artery  having  a  lumen  diameter  of  0.4  cm  and 
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Figure  21.  SAR  Distribution  for  Model  1 
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second  of  heat  generation. 


ire  23.  Two-dinen8ional  contour  plot  of  the  teaperature  distribution  for  Model  1  after 
second  of  heat  generation. 


Theme  1  Distribution  For  Model  2  (ot  t= 1  sec) 
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wall  thickness  of  1  am.  The  Maximal  thickness  of  the  plaque  layer  used 
in  this  Model  was  0.5  mm.  The  Microwave  angioplasty  Model  was  imposed 
on  a  200  by  80  cylindrical  grid  systen  with  a  uniforM  cell  size  of  1  mm 
in  both  x  and  r  directions.  To  Meet  stability  criteria,  the  therMal 
transient  solution  was  obtained  with  a  time  step  of  2.5  x  10*^  seconds. 
Although  the  smaller  cell  size  increased  the  number  of  iterations  and 
coaputational  time  it  allowed  for  a  more  detailed  analysis  of  theraal 
diffusion  within  the  blood  vessel  wall.  A  diagram  depicting  the  third 
Model  is  shown  in  Figure  27.  The  initial  SAR  calculated  by  the 
electromagnetic  FDTD  program  and  final  temperature  distribution  results 
are  shown  in  Figures  28,  29,  and  30,  respectively. 


Figure  27.  Diagram  of  Model  3 
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rigure  29.  Three-disensional  surface  plot  of  the  tesperature  distribution  for  Model  3 
if ter  one  second  of  heat  generation. 


Therwnl  Distribution  For  Model  3  (at  t=1  sec) 


Figure  30.  Two-dimensional  contour  plot  of  the  temperature  distribution  for  Model  3  after 
one  second  of  heat  generation. 
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Analysis  of  Thermal  Results 

The  success  of  the  microwave  angioplasty  techniques  is  dependent 
on  three  Major  criteria.  The  first  criterion  is  to  heat  the 
atheromatous  plaque  to  a  high  enough  temperature  for  removal  by  Belting 
or  vaporization.  Most  biological  tissues  contain  high  concentrations 
of  water  relative  to  other  cheaicals.  Previous  studies  with  laser 
angioplasty  techniques  have  shown  that  tissue  ablations  starts  at 
temperatures  higher  that  100*C,  which  is  closely  associated  with  the 
boiling  point  of  water  (Welch  et  al.,  1985).  Tissue  ablation  occurs 
after  the  water  in  tissues  is  initially  vaporized  and  is  accompanied  by 
a  rapid  rise  in  the  plaque  temperature.  The  rapid  temperature  rise  in 
tissues  is  caused  by  a  reduction  in  the  specific  heat  and  conductivity 
as  the  tissue  is  dehydrated.  As  mentioned  previously,  there  are  three 
different  types  of  atheromatous  plaques.  Depending  on  the  type  of 
plaque  predominantly  present,  effective  ablation  in  tissue  occurs  in  a 
temperature  ranges  of  160  to  310*C  (Welch  et  al.,  1985). 

Second,  the  prevention  of  significant  thermal  injury  to  normal 
vascular  and  muscle  tissue  is  of  primary  concern  in  applying  the 
microwave  angioplasty  technique.  Thermal  tissue  damage  is  defined  as 
loss  of  biological  function  in  cells  or  extracellular  fluid  (Welch, 
1985).  Tissue  damage  is  dependent  on  the  amount  of  time  the  tissues 
are  exposed  to  increased  temperature  in  addition  to  their  temperature 
rise.  In  the  case  of  hyperthermia,  the  body  can  withstand  temperatures 
up  to  a  range  of  42  to  43 *C  for  a  period  of  60  minutes  before  thermal 
injury  occurs.  At  higher  temperatures  the  body’s  resistance  for 
thermal  damage  is  decreased  and  cellular  necrosis  occurs  at  much 
shorter  periods  of  exposure  to  abnormally  high  temperatures.  The 
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Mathematical  relationship  to  predict  the  thermal  damage  to  tissue  as  a 
time-dependent  rate  process  was  derived  by  Henriques  and  Moritz 
(Henriques  and  Moritz,  1947;  Welch,  1985). 


<*n(r»zjt?  -  A  exD 
dt  A  exp 


[ 


RT ( r , z , t ) 


] 


(3-1) 


where:  fi  =  daaage  function 
A  -  a  constant 

E  =  activation  energy  for  a  reaction 
R  =  universal  gas  constant 
T  =  absolute  temperature 
t  =  time 

By  numerically  integrating  this  equation,  the  thermal  damage  can 
be  predicted  with  experimental  evaluation  of  the  constant,  A,  and 
activation  energy,  E.  Although  these  values  have  not  been  determined 
for  the  microwave  angioplasty  model,  the  results  for  cellular  necrosis 
due  to  thermal  radiation  in  tissues  (Figure  31)  can  be  used  to 
determine  whether  thermal  injury  will  occur  in  blood  vessel  and  muscle 
tissue.  For  example,  skin  tissues  at  75*C  will  reach  the  threshold 
damage  of  1.0  within  10“^  second. 

Blood  temperature  increase  caused  by  the  absorption  of  thermal 
energy  induced  by  the  electromagnetic  field  also  deserves  some 
consideration.  Previous  studies  have  demonstrated  that  prolonged 
temperatures  greater  than  48 'C  will  result  in  hemolysis  of  the  blood 
and  increased  red  cell  osmotic  fragility  (Wilson,  Knauf,  and  Iserson, 
1987).  Furthermore,  the  body  is  an  excellent  thermoregulator  at 
temperatures  below  45 *C,  but  it  becomes  less  efficient  for  significant 
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Figure  31.  Rate  of  accumulation  of  damage  in  skin  tissue  as  a  function 
of  increase  in  tissue  temperature.  This  graph  is  used  to  estimate  the 
number  of  seconds  required  before  threshold  damage  occurs.  (Adapted 
from  Welch,  1985) 


increases  in  blood  temperature.  Since  the  amount  of  microwave  exposure 
is  anticipated  to  be  extremely  short  fraction  of  a  second,  hemolysis  of 
the  blood  may  not  be  detrimental  to  the  welfare  of  the  patient.  On  the 
other  hand,  extremely  high  heat  generation  in  blood  may  cause  a 
substantial  imbalance  in  the  body’s  thermal  regulatory  system  if  the 
overall  body  temperature  is  raised  above  45 'C.  Therefore,  some 
evaluation  of  the  body’s  natural  ability  to  dissipate  the  generated 
heat  and  maintain  close  to  normal  body  temperatures  should  also  be 
considered. 

The  electromagnetic  fields  produced  by  the  microwave  applicator 
generated  one  to  three  sharp  spikes  of  extremely  high  SAR.  The  sharp 
increases  in  the  thermal  energy  absorbed  by  the  tissues  are  located 
either  near  the  area  surrounding  the  gap  in  the  outer  conductor  of  the 
microwave  applicator  or  where  the  electric  vector  field  intersects  the 
balloon  surface  in  a  normal  direction.  In  the  latter  case,  the 
difference  in  the  complex  permittivity  and  electrical  conductivity 
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properties  between  blood  and  air  results  in  a  significant  increase  in 
the  SAR  for  blood.  (A  more  thorough  discussion  of  this  concept  is 
covered  in  Appendix  A.)  Therefore,  the  occurrence  of  a  spike  is 
dependent  upon  the  aicrowave  applicator  design  and  the  shape  of  the 
electric  field.  The  location  of  this  spike  is  typically  at  the  ends  of 
the  balloon. 

The  temperature  distribution  after  heating  the  angioplasty  model 
for  one  second  yielded  a  heating  pattern  similar  to  the  SAR  output. 

Results  of  the  finite  difference  thermal  analysis  generally  showed 

atheromatous  plaque  heating  due  to  two  effects:  (1)  selective  heating 
and  (2)  convective  heat  transfer.  In  the  former  case,  plaque  was 
selectively  heated  by  the  microwave  applicator  in  the  area  closest  to 
the  gap  of  the  outer  conductor.  Convective  heat  transport  is  also 
present  because  of  the  fluid  flow  and  temperature  difference  between 
the  blood  and  plaque. 

In  Model  1,  plaque  heating  is  accomplished  by  three  SAR  spikes 
produced  by  the  microwave  applicator.  Consequently  three  areas  of  high 
heat  intensity  were  created.  In  each  of  these  three  areas,  most  of  the 

heating  occurred  in  the  blood  and  the  plaque  layers  closest  to  the 

blood.  One  second  of  heat  generation  in  the  plaque  resulted  in  peak 
plaque  temperatures  in  the  100 'C  range  for  the  SAR  spike  closest  to  the 
incoming  blood  flow,  200‘C  in  the  center,  and  65*C  at  the  downstream 
end  of  the  applicator.  In  all  three  cases,  blood  temperatures  were,  on 
the  average,  approximately  two  times  higher  than  the  plaque 
temperatures.  The  dissipation  of  heat  into  the  neighboring  blood 
vessel  resulted  in  relatively  high  vascular  temperatures  of  50 'C.  At 
the  downstream  end  of  the  balloon,  some  plaque  heating  is  caused  by 
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convective  heat  transfer  rather  than  selective  heating.  In  this  case, 
the  SAR  spike  generates  extremely  high  blood  temperatures  in  the  300  to 
400 ’C  range.  Although  the  convective  heat  transfer  raised  the  plaque 
temperatures  by  approximately  28 *C,  it  was  not  sufficient  to  promote 
ablation.  Figures  32  to  34  depicts  the  temperature  distributions  in 
each  of  the  three  SAR  spikes  as  a  function  of  the  radius,  r.  The 
temperatures  are  shown  for  the  blood,  plaque,  vascular,  and  muscle 
tissues,  excluding  the  microwave  applicator. 


Figure  32.  Temperature  distribution  for  various  materials  in  Model  1 
at  Peak  1.  Heat  generation  on  the  upstream  end  of  the  microwave 
applicator  was  caused  by  the  permittivity  difference  between  air  and 
blood.  Results  are  shown  for  a  one  second  microwave  pulse  at  t  =  1.0 
second. 
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Figure  33.  Temperature  distribution  for  various  materials  in  Model  1 
at  Peak  2.  Heat  generation  was  produced  by  the  microwave  applicator 
gap  located  in  the  middle  of  the  thermal  model.  Results  are  shown  for 
a  one  second  microwave  pulse  at  t  *  1.0  second. 


0.00929  0.01643  0.02358  0.03072  0.03787 

Radius,  m 

Figure  34.  Temperature  distribution  for  various  materials  in  Model  1 
at  Peak  3.  Heat  generation  on  the  downstream  end  of  the  microwave 
applicator  was  caused  by  the  permittivity  difference  between  air  and 
blood.  Results  are  shown  for  a  one  second  microwave  pulse  at  t  *  1.0 
second. 
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Angioplasty  Model  2  produced  a  single  SAR  spike.  Unlike  the 
first  model,  temperatures  in  the  atheromatous  plague  after  one  second 
of  heat  generation  were  not  high  enough  to  initiate  ablation.  Most  of 
the  thermal  energy  was  absorbed  by  the  blood  rather  than  plaque.  Blood 
temperatures  ranged  from  300  to  380 ‘C  with  corresponding  plaque 
temperatures  of  approximately  70  to  97 ‘C  in  the  area  of  the  SAR  spike 
location.  Despite  the  lack  of  adequate  heating  in  the  plaque,  blood 
vessel  temperatures  also  increased  to  about  60  to  67 *C.  The 
temperature  distribution  at  the  site  of  peak  SAR  for  Model  2  is  shown 
in  Figure  35. 

Although  two  SAR  spikes  were  generated  by  the  microwave 
applicator  in  the  third  model,  only  one  provided  any  significant 
heating  in  either  the  blood  or  atheromatous  plaque.  Much  of  the  heat 
generated  in  the  blood  by  the  first  spike  is  cooled  by  the  incoming 
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Figure  35.  Temperature  distribution  for  various  materials  at  peak  heat 
generation  for  Model  2.  Results  are  shown  for  a  one  second  microwave 
pulse  at  t  =  1.0  second. 
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blood  flow.  Blood  temperature  increases  were  kept  to  a  minimum  as 
indicated  by  relatively  low  average  blood  temperatures  of  40  to  45 *C 
for  this  region.  Some  thermal  energy  was  also  absorbed  in  the  blood 
vessel  and  muscle.  Therefore,  undesirable  temperature  increases  in  the 
vascular  tissues  were  seen  next  to  the  entrance  region.  The  average 
tissue  temperatures  in  this  region  were  in  the  70 *C  range. 

The  second  SAR  peak  in  Model  3  provided  sufficient  thermal  heat 
generation  in  the  atheromatous  plaque  to  cause  ablation  to  occur.  The 
microwave  applicator  was  successful  in  delivering  more  heat  generation 
in  the  plaque  layers  instead  of  the  blood.  The  maximum  plaque 
temperature  in  this  region  of  the  model  was  181 *C.  Temperatures  in  the 
blood  were  lower  than  those  for  plaque  by  an  average  of  50 *C.  Although 
this  model  was  more  successful  in  achieving  selective  plaque  ablation, 
extremely  high  tissue  temperatures  were  found  in  the  vascular  tissue. 
The  amount  of  heat  in  the  1.0  mm  thick  blood  vessel  was  conducted  away 
quite  rapidly.  For  example,  the  temperature  of  the  blood  vessel  cell 
next  to  the  plaque  was  127 *C.  For  the  outermost  blood  vessel  cell,  the 
temperature  dropped  to  46.1*C.  A  temperature  distribution  for  this 
region  of  the  model  is  shown  in  Figure  36. 

Because  many  simplifying  assumptions  were  made  to  model  the  fluid 
flow  and  convective  heat  transfer  in  this  thermal  analysis,  a  simple 
numerical  experiment  was  conducted  to  determine  the  importance  of  this 
contribution.  The  subroutine  to  calculate  the  convective  heat  transfer 
coefficients  was  bypassed  and  all  coefficients  were  set  to  zero.  As 
seen  in  Figure  37,  the  blood  temperatures  for  the  case  with  convection 
nearly  overlaps  the  temperature  distribution  without  convective  heat 
transfer  after  one  second  of  constant  microwave  heat  generation.  There 
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Figure  36.  Temperature  distribution  of  various  materials  at  peak  heat 
generation  for  Model  3.  Results  are  shown  for  a  one  second  microwave 
pulse  at  t  =  1.0  second. 


Figure  37.  Comparison  of  the  effects  of  convective  heat  transfer  in 
blood  temperatures  using  Model  1.  Results  are  shown  for  a  one  second 
pulse  at  t  =  1  second.  Both  curves  overlap  in  most  areas  of  the  curve. 
Therefore,  convection  is  not  of  significant  importance  during  short 
periods  of  microwave  heating. 


is  a  ssall  region  in  the  center  where  a  small  difference  exists, 
however,  this  variation  did  not  alter  the  results  of  this  thermal 
analysis  during  short  time  periods  of  constant  heat  generation. 

As  Model  1  was  allowed  to  cool  for  three  additional  seconds  after 
one  second  of  heat  generation,  an  increased  difference  was  observed 
between  the  case  with  convective  heat  transfer  and  the  same  model 
calculated  without  convection.  These  results  are  graphically  shown  in 
Figure  38.  These  results  demonstrate  the  importance  of  developing  a 
better  numerical  aodel  for  the  blood  flow,  especially  during  the 
cooling  process  and  possibly  between  microwave  pulses. 

In  all  of  the  models  analyzed,  blood  and  vascular  tissue 
temperatures  rose  above  the  limits  for  preventing  thermal  injury.  As  a 
possible  resolution  to  this  problem,  a  limited  study  of  microwave 
pulsing  was  also  conducted.  By  inducing  a  short  high  microwave  energy 
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Figure  38.  Comparison  of  the  effects  of  convective  heat  transfer  in 
blood  temperatures  using  Model  1  after  a  one  second  pulse  of  heat 
generation  plus  three  additional  seconds  of  cooling  ( t  =  4  sec).  After 
a  few  seconds  of  cooling  the  convective  heat  transfer  due  to  fluid  flow 
becomes  more  significant  in  the  thermal  analysis. 
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input  to  the  system  and  turning  it  off,  temperatures  high  enough  to 
initiate  ablation  could  be  achieved  while  letting  the  thermal  system 
cool  to  avoid  permanent  theraal  damage. 

To  test  the  thermal  effect  of  aicrowave  pulsing,  a  0.20  second 
pulse  was  applied  to  Model  3.  The  intensity  of  the  SAR  distribution 
was  doubled  in  order  to  achieve  plague  ablation  within  0.20  seconds. 

The  aodel  was  allowed  to  cool  up  to  3  seconds.  The  result  of  this 
analysis  is  exhibited  in  Figures  39  and  40.  The  temperatures  in  the 
blood  vessel  layers  increased  because  the  direction  of  heat  dissipation 
is  toward  the  cooler  muscle  tissue.  Since  temperatures  in  the  blood 
were  usually  higher  than  those  for  plaque,  more  heat  was  being  added  to 
the  plague  which  eventually  conducted  into  the  blood  vessel  and  muscle 
layers.  Temperatures  in  the  vascular  tissues  did  exceed  45 ‘C 
especially  in  the  layers  next  to  the  plague  and  blood  vessel  interface 
due  to  conduction  into  the  muscle  tissue.  These  temperatures  were 
slightly  higher  immediately  following  heat  generation.  For  example,  at 
a  specific  location  in  the  blood  vessel  having  a  temperature  of  60.4’C 
after  a  0.20  second  pulse  of  microwave  heat  generation,  the  temperature 
rose  up  to  68.5'C  one  second  after  turning  off  the  heat  source.  Hot 
spots  caused  by  the  electrical  field  crossing  the  air-blood  interface 
created  significant  heating  of  the  blood.  In  Model  1  blood 
temperatures  reached  as  high  as  245‘C.  Although  blood  temperatures  in 
Model  3  at  the  peak  site  of  heat  generation  were  less  than  plaque 
temperatures,  the  amount  of  heat  convected  to  blood  was  minimal  in 
comparison  to  conduction  in  the  tissues. 

The  thermal  properties  of  atheromatous  plaques  used  in  all  models 
discussed  above  were  based  on  a  combined  average  of  property  value  for 
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Figure  39.  Heat  dissipation  after  a  0.20  second  pulse  of  microwave 
heat  generation  using  Model  3.  Temperature  distributions  are  shown  for 
ts0.20  sec,  t=l  sec,  and  t=3  sec. 
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Figure  40.  Temperature  distribution  in  tissues  after  a  0.20  second 
pulse  of  microwave  heat  generation  using  Model  3.  Temperatures  are 
shown  for  t=0.20  sec,  t=l  sec  and  t=3  sec  as  in  Figure  39. 


fatty,  fibrous,  and  calcified  plaques  (Welch,  1985).  In  different 
parts  of  the  body,  certain  predominant  types  of  plaques  are  more  likely 
to  be  deposited  on  the  blood  vessel  linings.  For  example,  calcified 
plaques  are  more  commonly  found  in  the  blood  vessels  of  the  leg  and 
limbs  whereas  fatty  and  fibrofatty  plaques  are  usually  deposited  in  the 
coronary  arteries.  Further  numerical  experiments  were  performed  to 
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determine  whether  certain  types  of  plaques  would  have  a  higher  affinity 
for  selective  plaque  ablation. 

The  amount  of  thermal  energy  absorbed  in  the  plaque  varied 
according  to  the  predominant  plaque  type  in  the  angioplasty  model. 
Because  of  the  higher  water  content,  the  thermal  conductivity  and 
specific  heat  values  for  fatty  plaques  were  the  highest.  Applying  an 
equal  amounts  of  microwave  heat  generation  to  Model  3  for  all  three 
cases  resulted  in  consistently  higher  temperature  rises  for  fatty 
plaques.  As  shown  in  Figure  41  fatty  plaques  are  followed  by  fibrous 
then  calcified  types  in  their  ability  to  absorb  thermal  energy.  The 
variance  between  plaque  types  is  not  enormous.  At  site  of  peak  plaque 
temperature  the  fatty  plaque  was  about  15.5  percent  higher  than  the 
calcified  type.  The  difference  was  much  less  for  fibrous  plaque  with  a 
three  percent  difference  in  comparison  to  the  fatty  plaque.  Therefore, 
application  of  this  angioplasty  technique  is  not  limited  to  the  removal 
of  fatty  plaques  alone.  Plaques  normally  consist  of  a  mixture  of 
various  plaque  types.  Since  the  variance  in  results  was  not 
significant,  the  thermal  analysis  averaging  all  three  types  of  plaques 
was  not  a  bad  approximation  for  calculating  the  overall  plaque  thermal 
properties. 
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Figure  41.  Temperature  distributions  for  various  plaque  types.  Fatty 
plaques  absorb  less  than  15.5  percent  more  thermal  energy  than  other 
types.  Results  are  shown  for  a  one  second  microwave  pulse  at  t=l  sec. 


Although  the  results  of  this  thermal  analysis  have  not  been 
checked  experimentally,  similar  results  were  obtained  in  ablation 
studies  for  the  laser  angioplasty  technique.  In  research  work  done  by 
Welch  et  al.  a  laser  probe  was  applied  to  an  atherosclerotic  human 
aorta  with  a  wall  thickness  of  1.5  mm.  The  atheromatous  plaque  was 
heated  up  to  180 ‘C  to  induce  ablation  and  the  adventitial  temperatures 
reached  70*C  during  the  ablation  process  (Welch  et  al.,  1987).  The 
heat  transport  mechanisms  for  both  laser  and  microwave  angioplasty 
techniques  are  comparable  since  both  methods  attempt  to  heat  the  plaque 
to  the  same  temperature  levels.  Heat  in  both  cases  is  dissipated 
radially  outward  toward  the  cooler  muscle  tissue.  The  blood  vessel 
temperatures  in  the  thermal  models  analyzed  by  this  study  reached  the 
60  to  80‘C  range.  In  Model  3  a  maximum  temperature  of  127 'C  was 
observed  at  the  plaque  and  vessel  wall  interface.  The  similarities  in 
these  analyses  provides  a  fair  amount  of  confidence  in  the  validity  of 
this  thermal  model. 

During  the  development  of  the  finite  difference  code  several 
adjustments  and  approximations  were  made  to  allow  the  use  of  SAR  data 
generated  by  the  Department  of  Electrical  Engineering.  For  one,  the 
data  obtained  provided  relative  SAR  data  rather  than  actual  heat 
generation  values.  Results  from  the  FDTD  program  produced  SARs  having 
an  order  of  magnitude  of  10“®  for  heat  generated  by  the  microwave. 

Over  a  wide  range  of  the  microwave  exposures,  the  SARs  linearly 
increase  with  increasing  power  input.  Therefore,  to  show  effective 
heating  in  the  thermal  model,  a  multiplier  was  estimated  to  increase 
SAR  values  to  a  level  where  the  plaque  ablation  could  take  place. 


54 

Further  studies  to  include  the  plaque  ablation  portion  of  finite 
difference  code  will  require  a  correction  of  this  discrepancy  to  allow 
simultaneous  execution  of  both  FDTD  and  FDATA  programs. 

Additionally,  the  relative  SAR  data  ware  calculated  at  the  center 
of  each  finite  diiference  cell  rather  than  the  boundary  interfaces. 
Because  both  conductive  and  convective  phenomena  takes  place  in  this 
thermal  analysis,  the  SAR  data  required  recalculations  to  shift  the 
coordinate  system  by  a  half  interval  in  both  the  x  and  r  directions. 
Therefore,  a  cubic  spline  interpolation  subroutine  the  SAR  values  was 
implemented  to  calculate  this  change.  Models  1  and  2  were  originally 
generated  by  the  Department  of  Electrical  Engineering  to  test  various 
applicator  designs  numerically.  The  coordinate  mesh  in  these  to  models 
was  too  coarse  (0.1429  cm)  and  the  cubic  spline  interpolation  produced 
unrealistic  numbers  where  large  fluctuations  in  the  SAR  data  were 
observed.  This  problem  was  eliminated  by  using  a  linear  interpolation 
subroutine  instead.  Therefore,  a  5  to  10  percent  error  in  this 
solution  may  be  present  as  a  result  of  the  coordinate  system  changes 
between  the  FDTD  and  FDATA  programs. 

It  is  important  to  note  that  Models  1  and  2  do  not  reflect  the 
actual  dimensions  of  the  vascular  system.  For  example,  the  model  data 
shows  the  blood  vessel  starting  15  cells  away  from  the  centerline.  For 
the  cell  size  of  0.1429  cm  in  the  r  direction  this  equates  to  a  4  cm 
lumen  diameter.  The  largest  blood  vessel  in  the  human  body  is  2.5  cm 
in  diameter.  The  blood  vessel  was  I'lo  modelled  as  being  only  one  cell 
thick.  Model  3  was  developed  to  account  for  this  discrepancy  and  to 
better  evaluate  the  heat  distribution  within  the  vascular  tissues. 
Therefore,  Models  1  and  2  are  useful  in  determining  whether  or  not 
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selective  plague  ablation  can  be  achieved  but  they  do  not  represent  the 
actual  physical  system. 

In  Model  3  the  BAR  spike  at  the  entrance  region  of  the  model  was 
not  anticipated.  The  insulative  boundary  conditions  in  this  case  did 
not  accurately  represent  semiinfinite  conditions  because  of  the  build 
up  of  heat  in  this  region.  Therefore,  the  high  adventitial 
temperatures  observed  in  this  study  are  probably  lower  than  the 
calculated  values.  Boundary  condition  problems  of  this  nature  can  be 
alleviated  by  changing  the  dimensions  of  the  cell  size.  By  generating 
a  cylindrical  grid  where  the  dx  dimension  is  two  or  three  times  greater 
than  dr,  the  semi infinite  boundary  conditions  can  be  met  without 
increasing  the  matrix  size  and  computational  iterations.  As  an 
alternative  method  for  modelling  the  boundary  conditions,  the  boundary 
temperatures  could  be  held  constant  at  normal  body  temperature  (37 ’C). 

The  extent  of  permanent  thermal  damage  to  the  vascular  tissues 
can  not  be  clearly  defined  at  this  time.  Other  thermographic  and 
histological  studies  with  laser  angioplasty  techniques  have 
demonstrated  that  actual  ablation  of  normal  vascular  tissue  requires 
temperatures  exceeding  180 'C.  Tissue  dehydration  and  shrinkage  starts 
at  temperatures  of  60*C  or  greater  (Welch  et  al. ,  1987).  For  small 
increments  of  time,  the  tissues  are  able  to  tolerate  different  levels 
of  heating  without  causing  irreversible  permanent  damage.  Therefore, 
although  high  temperatures  in  the  blood  vessel  tissues  were  observed  in 
this  thermal  analysis,  short  exposures  to  the  elevated  temperatures  may 
not  cause  a  significant  amount  of  thermal  destruction. 

In  this  thermal  analysis,  hemolysis  of  the  blood  appeared  to  be  a 
potential  problem.  However,  the  extent  of  adverse  thermoregulation 
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alteration  and  blood  loss  is  expected  to  be  minimal  for  short  exposures 
of  aicrowave  heat  generation.  The  circulatory  system  stores  an  average 
of  8  liters  of  blood  within  its  vast  system  of  arteries,  veins  and 
capillaries.  The  aaount  of  heated  blood  is  anticipated  to  be  far  less 
than  the  total  blood  voluae  available.  This  assuaption  holds  true  for 
the  laser  angioplasty  technique.  With  iaproveaents  to  the  present 
applicator  design  to  eliminate  the  hot  spots  in  the  blood  created  by 
the  air-blood  interface,  blood  hemolysis  should  not  be  a  major  concern. 

Pulsing  applications  of  the  microwave  angioplasty  technique 
provided  better  results  than  the  constant  heat  generation  models.  High 
plaque  heating  for  ablation  could  be  accomplished  while  minimizing 
temperature  rises  in  the  blood  vessel  and  muscle  tissues.  Further 
numerical  optimizing  experiments  with  high  heat  generation  at  multiple 
time  exposures  may  provide  more  desirable  temperature  profiles. 

This  thermal  analysis  was  a  first  attempt  for  approximating  the 
conductive  and  convective  heat  transfer  contributions  to  the  microwave 
angioplasty  model.  Increasing  the  complexity  of  this  analysis  by 
eliminating  many  of  the  simplifying  assumptions  will  greatly  improve 
the  accuracy  of  the  results  obtained  here. 

As  previously  mentioned,  tissue  perfusion  and  metabolic  heat 
generation  were  omitted  from  this  thermal  analysis.  In  the  muscle 
tissue  the  omission  of  tissue  perfusion  may  not  be  such  a  good 
approximation  since  this  is  the  body's  only  means  of  heat  removal.  An 
inclusion  of  the  perfusion  term  in  the  bioheat  equation  will  allow  the 
heat  to  dissipate  at  a  faster  rate  than  observed  for  the  blood  vessels 
of  this  first  approximation. 
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A  better  treatment  of  the  fluid  flow  and  solution  for  the 
convective  heat  transfer  coefficient  is  also  required.  Preliminary 
observations  showed  some  significance  of  convective  heat  transfer  over 
long  periods  of  time  during  cooling.  These  results  should  be  verified 
by  conducting  a  more  detailed  analysis  of  the  blood  flow.  Further 
development  of  the  thermal  model  should  include  non-Newtonian  and 
pulsating  blood  flow  effects  as  well  as  account  for  the  curvature  and 
elasticity  in  the  blood  vessels.  Because  of  the  complex  geometry  of 
the  system,  the  fluid  flows  do  not  have  fully  developed  velocity  and 
temperature  profiles.  Entrance  region  effects  should  be  carefully 
examined  and  included  in  the  model.  Development  of  a  numerical 
solution  to  take  these  effects  into  consideration  may  alter  the 
convective  heat  transfer  coefficients  used  in  the  analysis. 

In  addition,  plaque  deposition  in  arteries  is  generally 
nonuniform.  Therefore,  a  three-dimensional  thermal  analysis  is  also 
required  to  handle  the  eccentricity  of  plaque  distribution  within  the 
blood  vessel. 

Thermal  properties  in  tissues  are  not  constant  at  different 
temperatures  and  pressures  as  modelled  in  this  analysis.  Additional 
studies  should  make  some  provisions  for  estimating  the  actual  thermal 
properties  as  a  function  of  temperature.  This  is  highly  recommended 
for  the  blood  and  atheromatous  plaque  layers  where  the  greatest 
temperature  gradients  occur  in  the  model. 

A  phase  change  portion  of  the  thermal  analysis  should  be  added  to 
the  numerical  solution.  The  results  presented  here  provide  temperature 
distributions  with  the  ablated  plaque  still  in  place.  By  removing  the 
plaque  from  the  model  and  recalculating  the  SAR  distribution 
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accordingly,  temperatures  in  the  blood  vessels  may  be  lower  than  the 
result  obtained  in  this  study.  As  ablation  takes  place  the  latent  heat 
is  carried  away  by  the  blood  flow.  Therefore,  plague  temperatures 
would  be  cooler  when  the  heat  source  is  turned  off  and  less  heat  will 
be  conducted  to  the  blood  vessel.  Heat  will  be  convected  to  the  outer 
layers  of  plague  due  to  the  increased  blood  temperature,  but  this  heat 
input  is  expected  to  be  relatively  small. 

Extremely  high  temperatures  in  the  blood  were  not  initially 
anticipated.  However,  the  heat  intensities  observed  in  this  study 
indicated  that  blood  vaporization  would  probably  occur  in  a  real  model. 
Some  treatment  of  the  phase  change  in  blood  is  needed  to  increase  the 
accuracy  of  the  thermal  model. 

Additional  thermal  studies  are  reguired  for  improving  the 
modelling  technigues.  For  example,  the  ablation  properties  of 
atheromatous  plagues  are  not  known.  Future  experimental  work  should 
include  measurements  for  the  latent  heat  of  fusion,  vaporization,  and 
sublimation  in  atheromatous  plagues  and  is  essential  to  accurately 
model  the  ablation  aspect  of  the  thermal  model.  The  rate  at  which  the 
plague  melts  or  vaporizes  also  reguires  evaluation. 

In  addition,  further  studies  to  determine  the  thermal  destruction 
thresholds  for  vascular  tissues  are  needed  to  evaluate  the  feasibility 
for  applying  this  angioplasty  technigue  clinically.  Studies  similar  to 
those  accomplished  for  skin  burns  are  reguired  to  calculate  burn  injury 
criteria  and  to  find  a  correlation  for  the  damage  constant  as  a 
function  of  time. 

Although  the  results  in  this  s^udy  show  vascular  tissue  damage 
and  hemolysis  might  have  occurred  during  treatment,  further  development 
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in  applicator  design  say  alleviate  these  problems.  Design  improvements 
for  the  microwave  applicator  should  include  a  heat  exchanger  system  for 
the  blood.  By  keeping  blood  temperatures  lower  than  those  of  the 
plaque,  heat  flows  into  the  blood  vessel  will  be  reduced  if  heat  can  be 
convected  away  by  the  blood.  Injection  of  blood  cooler  than  normal 
body  temperature  during  microwave  heating  may  help  to  minimise  blood 
temperature  increases. 

The  microwave  applicator  design  should  also  include  a  means  for 
directing  the  electromagnetic  field  better.  Plaque  distributions  in 
blood  vessels  are  not  axisymmetric.  Future  microwave  applicator 
designs  should  be  capable  of  heating  the  blood  vessels  in  sections  to 
avoid  heating  of  the  blood  vessel.  Undesirable  heating  at  the  air- 
blood  interface  of  the  applicator  ends  should  be  eliminated  by 
optimizing  the  applicator  balloon  shape  with  the  solution  of  the 
electromagnetic  field. 

Of  the  many  deficiencies  identified  here,  the  most  important 
improvement  is  the  quantifying  of  ablation  characteristics  in  the 
atheromatous  plaques  through  in  vitro  experimental  studies.  The 
success  of  the  microwave  angioplasty  method  is  entirely  dependent  upon 
the  ablation  process  in  plaques.  Adequate  atheromatous  plaque  removal 
depends  on  the  rate  in  which  the  plaque  melts  or  vaporizes  as  well  as 
the  latent  heat  required  to  initiate  a  phase  change.  Therefore,  to 
accurately  simulate  the  heat  transfer  in  an  angioplasty  model,  the 
phase  change  phenomena  must  be  added  to  the  thermal  analysis.  In  this 
study  temperatures  were  allowed  to  rise  constantly  as  heat  dissipated 
into  the  vascular  and  muscle  tissues.  Differences  in  temperature 
distributions  for  studies  including  the  ablation  phenomena  are 
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anticipated  since  blood  and  plaque  temperatures  are  held  constant 
during  the  phase  change  process.  The  estimated  power  required  to 
remove  the  remaining  plaque  through  a  series  of  subsequent  microwave 
pulses  will  also  be  affected  by  ablation.  These  experimental  studies 
will  moreover  provide  a  means  of  validating  the  thermal  results 
obtained  numerically.  Further  development  of  the  microwave  applicator 
design  would  be  meaningless  without  additional  analyses  of  the  plaque 
ablation  process. 

The  future  of  the  microwave  angioplasty  technique  is  promising. 
Implementation  of  the  improvements  described  above  may  lead  to  the 
practical  application  of  this  technique  clinically.  Thermal 
conductivity  from  the  plaque  to  blood  vessel  may  present  a  problem  if 
the  blood  vessel  cells  are  found  to  be  less  resistant  to  biological 
dysfunction  due  to  abnormal  temperature  rises.  The  application  of  this 
numerical  thermal  analysis  will  assist  in  predicting  the  occurrence  of 
thermal  injury  for  new  applicator  designs.  In  the  clinic  it  may  be 
useful  in  helping  physicians  estimate  the  microwave  dosimetry  required 
to  achieve  safe  and  effective  atheromatous  plaque  ablation. 


CHAPTER  4 


CONCLUSION  AND  RECOMMENDATIONS 

Conclusion 

In  this  study  a  finite  difference  thermal  sodel  was  developed  as 
a  first  approximation  to  predict  the  temperature  distribution  in  three 
percutaneous  transluminal  microwave  angioplasty  models.  The 
temperature  distributions  in  each  model  were  examined  to  determine 
whether  successful  atheromatous  plague  ablation  could  be  achieved 
without  causing  extensive  thermal  injuries  to  surrounding  vascular 
tissue  and  muscle. 

The  results  of  this  thermal  analysis  show  that  the  microwave 
angioplasty  technique  has  a  good  potential  for  application  in  a 
clinical  environment  in  the  treatment  of  atherosclerosis.  Selective 
atheromatous  plaque  heating  was  achieved  initially  and  temperatures 
well  above  100‘C  were  attained  for  plaque  ablation. 

High  temperatures  in  the  surrounding  tissues  and  blood  indicated 
a  high  probability  for  the  onset  of  thermal  injuries  during  application 
of  the  microwave  angioplasty  technique.  Temperatures  in  the  vascular 
tissues  did  exceed  45 *C  especially  in  the  layers  next  to  the  plaque  and 
ulouu  vessel  interface  due  to  conduction  into  the  muscle  tissue.  These 
temperatures  were  slightly  higher  immediately  following  heat 
generation.  Vascular  tissue  temperatures  were  in  the  60  to  80 ‘C  range 
one  second  after  the  microwave  applicator  was  turned  off.  Some  thermal 
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daaage  such  as  dehydration  and  shrinkage  to  the  tissues  would  occur  at 
this  elevated  temperature  level,  however,  temperatures  of  180 ‘C  or 
greater  are  needed  to  initiate  tissue  vaporization.  In  addition,  there 
is  a  high  potential  for  hemolysis  to  occur  in  the  blood  as  a  result  of 
the  hot  8 pots  caused  by  the  elw trical  field  crossing  the  air-blood 
interface  created  significant  heating  of  the  blood.  In  Model  1  blood 
temperatures  reached  as  high  as  245 *C. 

The  microwave  angioplasty  technique  performed  best  in  fatty 
plaques.  Higher  temperatures  were  observed  in  the  fatty  plaque  model 
for  identical  amounts  of  heat  generation  applied  to  all  three  plaque 
types.  However,  a  difference  in  maximum  temperatures  between  all  three 
plaque  types  was  not  prominent.  Therefore,  the  microwave  angioplasty 
technique  will  be  of  comparable  effectiveness  in  removing  all  plaque 
types. 


Recommendations  For  Further . Study 
Although  this  thermal  analysis  provides  a  good  first 
approximation  of  the  heat  transfer  phenomena  present  in  the  microwave 
angioplasty  system,  many  improvements  are  needed.  Further  studies  to 
analyze  the  phase  change  process  in  atheromatous  plaques  are  most 
important  since  ablation  is  the  primary  means  of  plaque  removal  using 
this  technique.  Additionally,  to  provide  a  more  exact  solution  of  the 
thermal  model  many  of  the  simplifying  assumptions  made  here  should  be 
eliminated.  The  results  previously  discussed  showed  that  convective 
heat  transfer  does  have  an  important  effect  on  the  thermal  analysis 
over  long  periods  of  cooling  and  potentially  between  microwave  pulses. 
A  more  accurate  method  of  analyzing  the  blood  flow  in  the  problem  is 
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therefore  required.  Further  studies  should  include  the  non-Newtonian, 
pulsating,  and  entrance  region  fluid  flow  effects.  The  solution  should 
also  account  for  curvature  and  elasticity  of  the  blood  vessels. 
Additional  studies  to  include  the  heat  transfer  effects  due  to  tissue 
perfusion  and  running  the  phase  change  version  of  the  FDATA  program 
will  greatly  improve  the  accuracy  of  these  results. 

Other  indirect  studies  are  needed  to  evaluate  the  safety  and 
possible  clinical  application  of  this  technique.  These  studies  should 
include  experimental  studies  of  phase  change  properties  in  plaque  and 
thermal  destruction  thresholds  for  vascular  tissues. 

Further  development  of  the  applicator  design  may  also  provide 
more  controlled  heating  in  the  plaque  layers  alone.  Examples  of  these 
improvements  consist  of  adding  a  heat  exchanger  system  to  regulate 
blood  temperatures,  better  applicator  design  to  prevent  undesirable 
heating  at  the  air-blood  boundary,  and  a  means  of  directing  the 
electromagnetic  field  for  non-axisymmetric  thermal  models. 


APPENDIX  A 


THE  MICROWAVE  ANGIOPLASTY  METHOD 


Microwave  Applicator  Design 

The  microwave  angioplasty  technique  for  treating  atherosclerosis 
is  siailar  to  the  balloon  and  laser  Methods  since  all  three  Methods 
eMploy  a  catheter  systeM  to  deliver  the  plaque  reaoval  device.  The 
Microwave  applicator  located  at  the  tip  of  the  catheter  consists  of  an 
extresely  fine  and  flexible  coaxial  cable  encased  in  a  thin  Metal  tube 
and  shorted  at  the  end  of  the  applicator.  The  coaxial  cable  is 
connected  to  an  external  Microwave  generator.  A  radial  gap  in  the 
Metal  casing  close  to  the  end  of  the  applicator  is  designed  to  create 
an  electrical  potential  difference  between  the  inner  coaxial  conductor 
(which  includes  the  coaxial  wire  and  shorted  end  of  the  tube)  and  the 
outer  coaxial  conductor  as  shown  in  Figure  42.  This  potential 
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Figure  42.  Microwave  applicator  as  designed  by  the  University  of  Utah 
Departsent  of  Electrical  Engineering. 
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difference  between  the  inner  and  outer  coaxial  conductors  produces  an 
electromagnetic  field.  The  surrounding  blood,  plague,  vascular,  and 
muscle  tissue  absorbs  the  thermal  energy  induced  by  the  high-energy 
electromagnetic  field  as  a  function  of  the  electrical  properties  of 
each  layer  of  tissue. 

Concept  of  Selective  Plaque  Ablation 
A  balloon  is  used  in  conjunction  with  the  microwave  applicator; 
however,  its  function  is  quite  different  than  its  use  in  the  coronary 
balloon  angioplasty  method  and  microwave  angioplasty  applicator 
proposed  by  Rosen  (1990).  In  hyperthermia  studies  conducted  by  the 
Christensen  and  Durney  (1981)  of  the  University  of  Utah,  it  was  found 
that  separating  the  fat-muscle  layer  and  electric  dipole  by  an  air  gap 
caused  the  fat  to  overheat.  Since  the  muscle  has  a  relatively  high 
electrical  conductivity  in  comparison  to  the  fat  and  muscle,  this 
causes  the  system  to  behave  as  two  capacitors  in  series  with  two 
dielectric  layers  as  shown  in  Figure  43. 


Figure  43.  The  electrical  conductivity  in  the  muscle-fat-air  system  is 
analogous  to  an  equivalent  electrical  circuit  with  two  capcitors  in 
series. 
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Assuming  the  electrical  and  thermal  properties  of  fat  are  similar  to 
that  of  atheromatous  plaques,  this  physical  phenomenon  is  applied  to 
the  microwave  angioplasty  model  to  selectively  heat  and  ablate  intra¬ 
vascular  atheromata  while  minimizing  thermal  injury  to  the  surrounding 
tissue. 

Selective  heating  of  the  fat  in  the  hyperthermia  model  also 
occurs  because  of  the  difference  in  the  complex  permittivity  of  the  fat 
and  muscle  tissues.  By  definition,  permittivity  is  the  electrical 
property  of  a  material  measuring  the  electromagnetic  field  interactions 
which  produce  thermal  energy  in  tissues.  These  interactions  include 
the  orientation  of  the  electric  dipoles  in  the  atoms  and  molecules  of 
the  tissues,  polarization  of  the  atoms  and  molecules,  and  displacement 
of  conductive  electrons  and  ions  in  tissues  (Christensen  and 
Durney, 1981 ) .  The  heat  generated  within  the  tissues  is  caused  by  the 
frictional  losses  due  to  the  movement  of  atoms  and  molecules  undergoing 
orientation  and  polarization  interactions  and  the  collisions  of  moving 
atoms  and  molecules  with  stationary  ones  during  the  displacement 
interaction. 

If  two  capacitor  plates  are  place  on  either  side  of  two  layers  of 
fat  and  muscle  as  shown  in  Figure  43,  a  strong  electric  field 
perpendicular  to  the  plates  is  produced.  A  much  weaker  magnetic  field 
is  also  produced  but  its  effect  is  not  significant  in  the  heat 
generation  of  the  tissues.  At  the  interface  of  the  fat  and  muscle,  a 
boundary  condition  (see  Figure  44)  relating  the  permittivities  of  each 
tissue  type  and  their  respective  electric  fields  is  given  as: 

e,E,  =  e.E. 


( A-l ) 
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Fat 

Muscle) 


Figure  44.  Boundary  condition  at  the  fat  and  muscle  tissues.  The 
plaque  having  similar  properties  of  fat  can  absorb  about  seven  times 
■ore  heat  generated  by  the  electromagnetic  field  than  the  muscle  and 
vascular  tissues. 


The  aaount  of  heat  absorbed  by  the  tissues,  also  known  as  the 
specific  absorption  rate  (SAR)  is  a  function  of  the  the  electrical 
conductivity  and  electromagnetic  field.  By  calculating  the  SAR  ratio 
of  fat  to  muscle  the  relationship  in  equation  (A-3)  is  obtained. 


SAR  = 


_  i  o!eL 


( A— 2 ) 
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By  substituting  the  boundary  condition  for  the  fat/auscle 
interface  into  equation  (A-3),  the  SAR  ratio  is  reduced  to  equation  (A- 
4).  Known  values  for  the  electrical  conductivities  and  coaplex 
permittivities  of  fat  and  muscle  are  substituted  into  this  equation 
giving  a  calculated  ratio  of  about  7. 


SAR, 

SAR. 


p.q 
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Published  values  for  the  electrical  properties  in  atheromatous 
plaques  are  unavailable.  However,  the  University  of  Utah  Hyperthermia 
Research  Group  was  able  to  obtain  experimental  measurements  for  plaque. 
Their  measurements  for  plaque  were  similar  to  that  of  documented  values 
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for  fatty  tissues.  For  exaaple,  the  experimental  relative 
permittivity,  e  and  electrical  conductivity,  a  values  for  plaque  were 
36  and  0.068  S/m,  respectively.  These  values  are  relatively  close  to 
the  values  for  relative  permittivity,  40  and  electrical  conductivity, 
0.019  S/m  in  fatty  tissues  (Hahn  et  al.,  1980).  Performing  the  same 
calculations  with  the  experimental  measurements  for  plaque,  the  SAR 
ratio  is  7.55  for  a  10  MHz  microwave  heat  source.  Both  documented 
electrical  properties  for  fatty  tissues  and  experimental  measurements 
for  plaque  demonstrates  the  high  potential  for  selective  heating 
because  fatty  tissues  absorb  approximately  seven  times  more  heat  than 
the  vascular  and  muscle  tissues. 


APPENDIX  B 


THERMAL  PROPERTIES  OF  BIOLOGICAL  TISSUES 


Tissue  Type 

Thermal 

Conductivity 

[W/(«‘C)] 

Average  Plaque  ( 1  ) 

0.4912 

Fatty  Plaque  (  1 ) 

0.4312 

Fibrous  Plaque  ( 1 ) 

0.4760 

Calcified  Plaque  (1) 

0.5667 

Blood  Vessel  Wall  (1) 

0.4320 

Muscle  (2) 

0.5100 

Blood  (3) 

0.6270 

Thermal 

Diffusivity  Heat  Capacity 

[■fysec]  [W*sec/(a^Bec) ] 


1. 265xl0“7 

3.878xl06 

1.241xl0~7 

3.476x10® 

1 . 295xl0~7 

3.677x10® 

1 . 364xl0-7 

4.480x10® 

1.220xl0*7 

3.541x10® 

1.800xl0“7 

2.833x10® 

1 . 545xlO~7 

4.059x10® 
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APPENDIX  C 


FINITE  DIFFERENCE  ANGIOPLASTY  THERMAL  ANALYSIS  PROGRAM 

Description  of  Nomenclature 

ALPHA  -  array  storing  the  thermal  diffusivities  for  each  tissue  type 
[a2/sec] 

C2  -  coaaon  coefficient  used  by  several  subroutine 
DR  -  differential  cell  radius  in  the  r  direction  [a] 

DT  -  differential  time  step  [sec] 

DX  -  differential  cell  length  in  the  x  direction  [a] 

HP  -  convective  heat  transfer  coefficient  for  the  plaque  [W/(a2‘C)] 

HM  -  convective  heat  transfer  coefficient  for  the  aicrowave  applicator 
[W/(a2,C)] 

I FIN  -  aaxiaua  nuaber  of  coluans  in  MODEL 

IOFF  -  integer  indicating  end  of  aicrowave  pulse  duration  INT(TIME/DX) 
IPRT  -  counter  for  interaediate  printouts  of  the  result 
ISTART  -  ainiaua  nuaber  of  coluans  in  MODEL 
JFIN  -  aaxiaua  number  of  rows  in  MODEL 

JRPL  -  integer  corresponding  to  the  location  where  the  plaque  starts 
for  a  given  x  location. 

JRPR  -  integer  corresponding  to  the  location  where  the  aicrowave  probe 
starts  for  a  given  x  location. 

JSTABT  -  ainiaua  nuaber  of  row  in  MODEL 

K  -  array  storing  theraal  conductivities  for  each  tissue  type  [W/a‘C] 
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MODEL  -  matrix  storing  the  model  parameters. 

1  -  refers  to  air 

2  -  refers  to  the  metal  of  the  microwave  applicator 

3  -  refers  to  the  muscle  tissue 

4  -  refers  to  the  atheromatous  plaque 

5  -  refers  to  the  blood 

6  -  refers  to  the  blood  vessel 

POWER  -  SAR  results  from  FDTD  program  [W/m3J 

QAVE  -  average  blood  flow  rate,  default  value  is  250  ml/min 

(91  -  matrix  storing  the  heat  generated  by  the  microwave  applicator 
[W/m3] 

QB  -  array  storing  the  lumped  heat  generated  in  the  blood  for  the  one 
dimension  convection  analysis  [W/m3] 

R0  -  radius  of  the  finite  difference  cell  closest  to  the  microwave 
applicator  [m] 

RHOC  -  array  storing  the  heat  capacity  (density  multiplied  by  the 
specific  heat)  for  various  tissue  types  [W»sec/(m3’C)] 

RPLX  -  array  storing  the  radii  values  at  the  plaque  layer  [m] 

RPRX  -  array  storing  the  radii  values  at  the  microwave  probe  [m] 

SEC  -  total  time  for  running  the  FDATA  program  [sec] 

SUM  -  summation  of  four  model  nodes  to  determine  which  finite 
difference  subroutine  to  execute 

TEMP  -  matrix  storing  the  initial  temperature  increases  computed  by  the 
FDTD  program  [ ’C) 

TEST  -  matrix  to  check  whether  the  program  is  executing  the  correct 
finite  difference  subroutine  for  the  plaque  geometry  given 

TIME  -  current  time  during  finite  difference  calculation 

TMAX  -  maximum  temperature  iound  in  the  final  results;  used  for 
normalizing  temperature  matrix  TN  and  allowing  its  use  with  EE  programs 

rc] 

TMP  -  matrix  storing  the  interpolated  temperatures  from  TEMP  [ *CJ 
TO  -  matrix  for  storing  the  old  temperature  distribution  [ ’C] 

TM  -  matrix  for  storing  the  new  temperature  distribution  [*CJ 

TP  -  array  for  storing  the  microwave  probe  surface  temperatures  [ *C] 

¥  -  average  blood  velocity  based  on  QAVE  [m/sec] 


1.  Input  the  total  number  of  seconds  to  evaluate  the  temperature 
distribution  in  the  angioplasty  model  (SEC). 

2.  Input  the  alternate  value  for  average  blood  flow  rate  in  ml/min  or 
use  the  default  value  250  ml/min. 

3.  Read  in  formatted  data  from  the  FDTD  program.  Data  includes  POWER, 
TEMP,  and  MODEL  matrices. 

4.  Initialize  all  temperature  matrices  and  arrays  (TO,  TB)  to  normal 
body  temperature,  37.0*C. 

5.  Interpolate  POWER  and  TEMP  to  compute  the  amount  of  heat  generation 
at  the  boundary  interfaces.  The  interpolated  data  are  stored  in  QM  and 
TMP  matrices. 

6.  Compute  the  radii  along  the  model  length  of  the  microwave 
applicator  and  plaque/blood  vessel  layer  using  MODEL.  Store  surface 
temperatures  for  the  nicrowave  applicator  form  TMP  into  TP. 

7.  For  each  finite  difference  cell  along  the  axial  direction,  compute 
the  amount  of  heat  in  the  blood  by  lumping  the  heat  generation,  QM  in 
the  radial  direction  per  unit  volume.  These  numbers  are  stored  in  QB. 

8.  Calculate  the  convective  heat  transfer  coefficients  using  cubic 
spline  interpolation  to  evaluate  points  between  values  given  in  Kays 
and  Crawford,  1980. 

9.  From  DX  to  SEC,  in  steps  of  DX,  do  the  following  steps. 

a.  Use  the  summation  of  MODEL(I,J),  MODEL(I+l,J),  MODEL(I , J+l ) , 
and  M0DEL(I+1,J+1)  to  determine  which  finite  subroutine  to  execute 
(SUM).  The  new  temperatures  calculated  from  matrix  TO  are  stored  in 
matrix  TN. 

i.  If  SUM=20,  run  R2MATL. 

ii.  If  SUM=18  or  22,  run  CONVBND  or  CONV3A/CONV3B. 

iii.  If  SUM= 1 6 ,  run  INTRIOR. 

iv.  If  SUM=17,  run  CONV5A/CONV5B. 

v.  If  SUM=19,  run  CONV4A/CONV4B. 

vi.  If  SUM=21,  run  CONV6A/CONV6B. 

b.  Compare  M0DEL(I,J)  with  MODEL  (I, J+l)  to  run  either  INTRIOR 
or  R2MATL  subroutines. 

c.  For  the  last  row  of  muscle  in  MODEL,  run  INSBND. 
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d.  Run  subroutines  CONVCNR,  INSCNR,  SI DEI ,  and  SIDE2  along  the 
edges  of  the  theraal  aodel. 

e.  Coapute  the  finite  difference  teaperatures  for  the  blood  flow 
using  subroutine  BFLOW. 

f.  Replace  TO  with  TN  for  the  next  iteration. 

10.  Noraalize  final  results  for  future  use  with  the  Departaent  of 
Electrical  Engineering  graphics  prograa  and  foraat  output. 

11.  Print  foraat ted  output. 


no  o  non  non 
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£******************************************************************** 

C  THESIS  PBOJECT  -  A  FINITE  DIFFERENCE  HEAT  TRANSFER  ANALYSIS  OF 
C  ATHEROMATOUS  PLAQUE  FOR  PERCUTANEOUS  TRANSLUMINAL  MICROWAVE 
C  ANGIOPLASTY  IN  BLOOD  VESSELS. 

C******************************************************************** 

FINITE  DIFFERENCE  ANGIOPLASTY  THERMAL  ANALYSIS  ( FDATA )  PROGRAM 

include  ’head.h’ 
parameter ( *=200 , n=200 , nt=6 ) 

REAL*8  TO(0:M,0: N) , ALPHA(6) ,RHOC(6) ,DX,DR,DT,TIME,QAVE, 

.  C2(6) ,K(6) ,TN(0:M,0:N) ,R0,COEFF,QM(0:M,0:N) ,TP(0:M) ,V(0:M) , 

.  TB(0:M+2) ,QB(0:M) ,RPLX(0:M) ,RPRX(0:M) ,HM(0:M) ,HP(0:M) ,SEC, 

POWER(M,N),TEMP(M,N),TMP(0:M,0:N),TMAX 
INTEGER  IPRT,MODEL(M,N) , JRPR(0:M) ,SUM, JFIN, IFN, IPRT2, IOFF, 

.  I TMP( M , N ) , I START ,  I FI N , J  START , TEST (0:M,0:N),JRPL(0:M) 

CHARACTER* 1  STATUS 
CHARACTER* 40  TITLE 

INITIALIZE  PARAMETERS 

COMMON  TO.TN.QM 

COMMON  /DAT/POWER , TEMP , MODEL , I START , I FI N , J ST ART , JFIN 
DATA  ALPHA ( 3 ) ,  ALPHA(6) ,  ALPHA( 4 ) ,  ALPHA (5  )/l .  8E-7 , 1 . 22E-7 , 

1 . 265E-7 , 1 . 54476E-7/ 

DATA  RHOC(3) ,RHOC(6) ,RHOC(4) ,RHOC(5)  /2.833E6.3.541E6, 

3. 878E6 , 4 . 05888E6/ 

DATA  K(3),K(6),K(4),K(5)  /. 51, .432, .4912,0.53/ 

DATA  DT , DX , DR  /2.5E-5.1.0000E-4.1.0000E-4/ 

IPRT=1 
IPRT2=1 
IFN=40 
IOFF=8000 

PRINT*, ’THIS  IS  THE  NONABLATION  VERSION  OF  THE  FDATA  PROGRAM.’ 
PRINT*,’  ’ 

INPUT  THE  NUMBER  OF  SECONDS  TO  EVALUATE  TEMPERATURE  CHANGE 

PRINT*, ’INPUT  THE  NUMBER  OF  SECONDS  TO  EVALUATE  TEMP  CHANGES’ 
READ*, SEC 
C 

PRINT*, ’DO  YOU  WANT  TO  INPUT  NEW  VALUES  FOR  BLOOD  FLOW  RATE?  TYPE 
.Y  TO  INPUT  NEW  VALUES  OR  N  TO  USE  A  DEFAULT  VALUE  OF  250  ML/MI N. ’ 
READ*, STATUS 
IF  (STATUS. EQ. ’Y’)  THEN 

PRINT*, ’INPUT  QAVE  IN  ML/MI N’ 

READ*,  QAVE 
ELSE 

QAVE=250.0 
END  IF 
C 

OPEN  ( 7 , FI LE= ’ ang9t . dat  * , status= ’ unknown ’ ) 

OPEN  ( 8 , FI LE=  *  aod9 . dat ’ , status= ’ unknown ’ ) 

OPEN  ( 9, FI LE=’ ang9fg.dat’ , status= ’ unknown ’ ) 
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OPEN  ( 10 , PI LE= ’ ang9 i . dat  * , status^ ’ unknown ’ ) 

C 

C  READ  IN  SPECIFIC  ABSORPTION  RATE,  TEMPERATURE,  AND  MODEL  DATA 
C  FROM  EE  FDTD  PROGRAM. 

C 

CALL  DATIN 

TITLES’ SPECIFIC  ABSORPTION  RATE  DATA’ 

CALL  PDATA2 ( 8 , POWER , TITLE ) 

TITLE= ’E-M  FIELD  TEMPERATURE  DATA  ’ 

CALL  PDATA2( 8, TEMP, TITLE) 

C 

C  ESTABLISH  INITIAL  TEMPERATURES  IN  THE  TO  MATRIX  AND  TB  ARRAY 
C  -  FOR  THIS  INITIAL  CASE,  ASSUME  ALL  TEMPERATURES  IN  EACH  LAYER  ARE 
C  THE  EQUAL  TO  THE  NORMAL  BODY  TEMPERATURE  OF  37  C 

C 

DO  10  I=ISTART-1 , IFIN 
TB( I )=37.0 

DO  20  JsJSTART-l.JFIN 
TO(I , J)=37.0 
20  CONTINUE 

10  CONTINUE 

TB( IFIN+1 )=TB( IFIN) 

TB(IFIN+2)=TB( IFIN+1 ) 

C 

DO  15  1=3, NT 

C2 ( I ) sCOEFF( ALPHA , DT , DX , I ) 

15  CONTINUE 

C 

C  SHIFT  POWER  MATRIX  BY  DR/2  AND  DX/2  TO  CALCULATE  THE  MICROWAVE  HEAT 
C  INPUT  AT  THE  BOUNDARY  AND  INTERFACE  NODES. 

C 

CALL  I NTRP1 ( 1 , QM , DR , DX ) 

TITLES 'INTERPOLATED  QM  DATA' 

CALL  PDATA1(10,QM, TITLE, TIME) 

CALL  I NTRP1 ( 2 , TMP , DR , DX ) 

TITLES ’INTERPOLATED  TEMP  DATA’ 

CALL  PDATA1(10, TMP, TITLE, TIME) 

CLOSE(UNIT=10) 

C  FIND  THE  RADII  ALONG  THE  MODEL  LENGTH  OF  THE  MICROWAVE  PROBE  AND 
C  PLAQUE  LAYER. 

C 

DO  30  I =1 START , IFIN 

DO  40  JsJSTART+1 , JFIN-1 

IF  ( (MODEL ( I , J ) . EQ. 5) . AND. ( MODEL ( I ,Jrl).NE.5))  THEN 
RPLX(I)=(J*DR) 

JRPL( I )sj 
ENDIF 

IF  ((MODEL(I ,J).EQ.5) .AND. (MODEL(I , J-l) .NE.5) )  THEN 
TP( I )=TMP( I ,  J-l ) 

JRPR( I )=J-1 
RPRX(I)sJRPR(I)*DR 
ENDIF 


noon  non  oooo  oooo  oooo 
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40  CONTINUE 
30  CONTINUE 

TP(0)=TP(1) 

ADJUST  PLAQUE  STARTING  J  POINT  TO  ACCOUNT  FOR  AN  INCREASE  IN  PLAQUE 
THICKNESS  AS  THE  BLOOD  APPROACHES  THE  BOUNDARY  LAYER. 

DO  41  I=ISTART, IFIN-1 

IF  (JRPL(I).GT.JRPL(I+1))  JRPL(I )=JRPL(I+1) 

41  CONTINUE 

FIND  THE  AVERAGE  RADII  IF  THE  COORDINATE  SYSTEM  IS  SHIFTED  TO  THE 
RIGHT  BY  DX/2. 

DO  35  I=ISTART, IFIN-1 

RPLX( I )=(RPLX( I )+RPLX(I*l ) )/2. 0 
RPRX( I )=(RPRX( I )+RPRX(I+l ) )/2.0 
35  CONTINUE 

RPLX ( 0 ) =RPLX ( 1 ) 

RPLX(IFIN)=RPLX( IFIN-1) 

RPRX ( 0 ) =RPRX ( 1 ) 

RPRX ( I  FI N ) =RPRX (IFIN-1) 

JRPL(0)=JRPL( 1 ) 

JRPR ( 0 ) = JRPR ( 1 ) 


LUMP  MICROWAVE  HEAT  INPUT  IN  THE  R  DIRECTION  FOR  EACH  LENGTH  OF  DX 
TO  REPRESENT  THE  TOTAL  HEAT  ADDED  TO  THE  BLOOD  PER  UNIT  VOLUME. 

DO  45  I =1 START, I FIN 

KK=( JRPL( I )-l )-( JRPR( I )+l ) 

IF  (KK.EQ.O)  THEN 
QB( I )=0. 5*QM( I , JRPL( I ) ) 

GO  TO  45 
ELSE 

DO  46  J=JRPR( I )+l , JRPL(I )-l 
Rl=DR*J+DR/2.0 
R0=DR*J-DR/2 . 0 

QB( I )=QB( I )+QM( I , J)*(R1**2-R0**2 ) 

46  CONTINUE 

R1=JRPL( I )*DR-DR/2 . 0 
R0=JRPR { I ) *DR+DR/2 . 0 

QB( I )=QB( I )/(Rl**2-R0**2 )+0. 5*QM( I , JRPR( I ) ) 

ENDIF 

45  CONTINUE 

QB(0)=QB( 1 ) 

COMPUTE  THE  AVERAGE  CONVECTIVE  HEAT  TRANSFER  COEFFICIENTS 
CALL  HTC( RPLX, RPRX, IFIN,HP,HM,K) 

DETERMINE  WHETHER  THE  NODE  BEING  EVALUATED  IS  AN  INTERIOR  NODE  OR 
ONE  WITH  BOUNDARY  CONDITIONS  AND  APPLY  THE  APPROPRIATE  SUBROUTINE. 


oooo  o  o  n  o 


77 


DO  25  TIM£=DT,SEC,DT 

TURN  OFF  HEAT  GENERATION  AT  THE  END  OF  MICROWAVE  PULSE  DURATION 

IF  (INT{TIME/DX) . EQ. IOFF)  THEN 
DO  48  11=0, IFIN 
QB(II}=0.0 
DO  47  JJ=0, JFIN 
QM(II ,JJ)=0.0 

47  CONTINUE 

48  CONTINUE 
ENDIF 

DO  50  I=ISTART,IFIN-1 
BO=JRPL(I)*DR 
DO  60  J=JRPR(I),JFIN-1 

IF  ( (MODEL(I , J) . EQ. 1) .OR. ( MODEL ( I , J ) .EQ. 2) )  GO  TO  60 

COMPUTE  THE  SUM  OF  THE  MODEL  NODES  TO  DETERMINE  WHICH  CONVECTIVE 
HEAT  TRANSFER  BOUNDARY  CONDITION  TO  APPLY  FOR  THE  PLAQUE. 

SUM=MODEL ( I , J ) +MODEL ( I , J + 1 ) +MODEL ( I + 1 , J ) +MODEL ( I + 1 ,  J + 1 ) 

IF  (SUM.EQ.20)  THEN 
IF(MODEL(I,J).EQ.4)  GO  Tn  62 
GO  TO  60 
ENDIF 

IF  ( { SUM.EQ. 18) .OR. (SUM.EQ. 22) )  THEN 
IF  ( MODEL ( I , J ) . EQ . 6 )  GO  TO  62 
IF  (MODEL(I , J) . EQ.MODEL( 1+1 , J ) )  THEN 
TEST ( I , J)=8 
IF  (SUM.EQ. 18)  TYPE=4 
IF  (SUM.EQ. 22)  TYPE=6 

CALL  CONVBND l HP , ALPHA , RHOC , DR , DT , TYPE , I , J , TB , C2 , RO ) 
R0=R0+DR/2.0 
GO  TO  60 
ELSE 

IF  (MODEL(I,J).EQ.5)  THEN 
TEST (I ,J)=9 

CALL  CONV3B ( HP , ALPHA , RHOC , DR , DX , DT , 4 , I , J , TB , C2 , RO ) 

GO  TO  61 
ELSE 

TEST( I , J)=10 

CALL  CONV3A ( HP , ALPHA , RHOC , DR , DX , DT , 4 , I , J , TB , C2 , RO ) 

GO  TO  61 
ENDIF 
ENDIF 
ENDIF 

IF  (SUM.EQ. 16)  THEN 
TEST( I , J )=1 

CALL  INTRIOR(ALPHA,RHOC,DR,DT,4,I,J,C2,RO) 

GO  TO  61 
ENDIF 

IF  (SUM.EQ. 17)  THEN 

IF  ( MODEL ( I , J ) . EQ . 5 )  THEN 


o  n  o  o 
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TEST(I,J)=11 

CALL  C0NV5B ( HP , ALPHA , RHOC , DR , DX , DT , 4 , I , J , TB , RO ) 

GO  TO  61 
ELSE 

TEST(I ,J)=12 

CALL  CONV5A(HP, ALPHA,RHOC,DR,DX,DT,4 , I , J ,TB,R0 ) 

GO  TO  61 
ENDIF 
ENDIF 

IF  (SUM.EQ.19)  THEN 

IF  (MODEL( I , J+l ) . EQ.  5 )  THEN 
TEST(I,J)=13 

CALL  CONV4A(HP,ALPHA,RHOC,DR,DX,DT,4,I,J,TB,C2,RO) 
RO=RO+DR/2 . 0 
GO  TO  60 
ELSE 

TEST( I , J)=14 

CALL  CONV4B(HP,ALPHA,RHOC,DR,DX,DT,4,I , J ,TB,C2,R0) 
R0=R0+DR/2 . 0 
GO  TO  60 
ENDIF 
ENDIF 

IF  (SUM.EQ.21)  THEN 

IF  ( MODEL ( I , J ) . EQ. 5 )  THEN 
TEST( I , J )=15 

CALL  CONV6B(HP,K,PHOC,DR,DX ,DT , 4 ,6, I , J ,TB,R0) 

GO  TO  61 
ELSE 

TEST( I , J )=16 

CALL  CONV6A(HP,K,RHOC,DR,DX ,DT,4 ,6, I , J ,TB,R0 ) 

GO  TO  61 
ENDIF 
ENDIF 

COMPUTE  NEW  TEMPERATURE  DISTRIBUTION  USING  THE  APPROPRIATE  THE 
INTERIOR  NODES. 

62  IF  (I.EQ.l)  THEN 

L=2 
ELSE 
L=I 
ENDIF 

IF  ( (MODEL( L-l , J ) . EQ.MODEL(L+l , J ) ) .AND. (MODEL(L, J) .EQ. 
MODEL(L, J+l ) ) )  THEN 
TEST( I , J )=1 

CALL  I NTRI OR ( ALPHA , RHOC , DR , DT , MODEL ( I , J ) , I , J , C2 , RO ) 

GO  TO  61 
ENDIF 

IF  ( MODEL (L,J ) .NE.MODEL( L, J+l ) . AND. (SUM. NE. 10) )  THEN 
TEST( I , J )=2 

CALL  R2MATL ( K , RHOC , DX , DR , DT , MODEL ( I , J ) , MODEL 
(I,J+1),I,J,R0) 

ENDIF 

61  R0=R0+DR 


noon 
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60  CONTINUE 
50  CONTINUE 

DO  70  I=ISTART,IFIN-1 
R0= JFI N*DR-DR/2 . 0 
TEST (I, JFIN) =3 

CALL  I NSBND ( ALPHA , RHOC , DR , DT , 3 , I , J , C2 , RO ) 

70  CONTINUE 
C 

C  COMPUTE  NEW  TEMPERATURE  DISTRIBUTION  AT  THE  ENDS  OF  THE  MODEL. 

C 

DO  90  I =1 START- 1 , 1  FI N , I  FI N 
IF  (I.EQ.O)  L=1 
IF  (I.EQ.IFIN)  L=IFIN 
R0=RPLX ( I ) 

DO  80  J=JRPL(I ) , JFIN 
IF  ( J . EQ. JFIN)  THEN 
TEST( I , J )=4 

CALL  INSCNR( ALPHA,RHOC,DR,DT , I , J , 3 ,C2 ,R0 , 1  FIN ,C2 ,R0 , 1  FIN ) 
GO  TO  80 
ENDIF 

IF  (MODEL(L, J ) . NE. 5 )  THEN 
IF  (MODEL(L,J).NE.MODEL(L,J+l))  THEN 
TEST ( I ,  J  )=5 

CALL  SIDE2(K,RHOC,DX,DR,DT,MODEL(L, J) ,MODEL(L, J+l) , 

.  I , J,R0, IFIN) 

GO  TO  81 
ENDIF 

IF  ( (MODEL(L, J ) . EQ.MODEL(L , J+l ) ) )  THEN 
TEST ( I , J)=6 

CALL  SIDE1 ( ALPHA, RHOC , MODEL ( L , J ) , DR , DT , , J ,C2 ,R0, IFIN) 

GO  TO  81 
ENDIF 
ELSE 

IF  (MODEL(L,J).EQ.MODEL(L,J+l))  GO  TO  80 
TEST( I , J )=7 

CALL  CON VCNR ( HP , ALPHA , RHOC , DR , DX , MODEL (L, J+l ) , I , 
J,TB,C2,R0,IFIN) 

R0=R0+DR/2.0 
GO  TO  80 
ENDIF 

81  R0=R0+DR 

80  CONTINUE 

90  CONTINUE 

IF  (TIME.EQ.DT)  THEN 
TITLE=’TEST  FILE’ 

CALL  ARR2(8 .TEST, M, 30, 100, 0,0, IFIN, JFIN, TITLE) 

CLOSE  (UNIT=8) 

ENDIF 

CALCULATE  THE  AVERAGE  BLOOD  FLOW  VELOCITY  THEN  CALCULATE  THE 
TEMPERATURE  DISTRIBUTION  OF  THE  BLOOD 

CALL  BFLOW ( TB , QB , RHOC , 5 , HP , HM , RPLX , RPRX , DX , DT , TP , V , 

.  IFIN,JRPR,JRPL,QAVE) 


o  o  o  o  o  o  n 
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STORE  NEW  TEMPERATURE  DISTRIBUTION  INTO  THE  OLD  TEMPERATURE  MATRIX 
FOR  THE  NEXT  ITERATION. 

DO  140  J=0, JFIN 
DO  130  1=0, I FIN 
TO( I , J)=TN( I , J ) 

130  CONTINUE 
140  CONTINUE 

PRINT  NEW  TEMPERATURE  DISTRIBUTION  FOR  THE  1000TH  TIME  STEP. 

IF  <IPRT.EQ.2000)  THEN 
I PRT= 1 

TITLES’ NEW  TEMP  PROFILE  ’ 

CALL  PDATA1 ( I FN,TN, TITLE, TIME) 

CLOSE  (UNIT=IFN) 

IFN=IFN+1 
ELSE 

IPRT=IPRT+1 
ENDIF 

25  CONTINUE 

C  NORMALIZE  FINAL  TEMPERATURE  DISTRIBUTION  OUTPUT  FOR  USE  WITH  THE  EE 
C  GRAPHICS  OUTPUT  PROGRAM.  (FIND  TMAX  FOR  TN(M,N) ;  DIVIDE  TN(M,N)  BY 
C  TMAX;  MULTIPLY  BY  A  FACTOR  OF  1000;  CONVERT  DATA  TO  AN  INTEGER 
C  MATRIX) 

C 

DO  350  I=ISTART-1 , IFIN-1 

nn  Qcn  1  —  TQT  APT— 1  TITTN-1 

TEMP( 1+1 , J+l )=(TN( I , J )+TN( 1+1 , J)+TN( I , J+l )+TN( 1+1 , J+l ) )/4.0 
360  CONTINUE 
350  CONTINUE 

TITLES’ NEW  TEMP  PROFILE  ’ 

CALL  PDATA2(7, TEMP, TITLE, TIME) 

CLOSE  (UNIT= 7) 

TMAX=0.0 

DO  310  I=ISTART-1,IFIN 
DO  320  JsJSTART-1 , JFIN 

IF  (TEMP( I , J ) .GT.TMAX)  TMAX=TEMP(I , J) 

320  CONTINUE 
310  CONTINUE 

DO  330  I =1 START, I  FIN 
DO  340  JsJSTART, JFIN 

ITMP(I , J )=INT( ( FLOAT ( TEMP ( I , J) )/FLOAT(TMAX ) )*1000) 

340  CONTINUE 
330  CONTINUE 

TITLE= ’FINAL  NORMALIZED  TEMPERATURE  OUTPUT’ 

CALL  ARRAYOUT (9,ITMP,M,30,100,I START , J  START , I FI N , JFIN, TITLE) 
TITLES’ MODEL  GENERATED  BY  EE’ 

CALL  ARRAYOUT ( 9 , MODEL , M , 30 , 1 00 , I START , JSTART ,IFIN, JFIN, TITLE) 
STOP 
END 


noon  ooooo 
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C  SUBROUTINE  TO  CALCULATE  T PRIME  (STORED  IN  MATRIX  X)  AT  DT  FOR  AN 
C  INTERIOR  NODE. 

C 

SUBROUT I NE  I NTR I OR ( AL PH A , RHOC , DR , DT , TYPE , I , J , C2 , RO ) 

INCLUDE  ’head. h’ 

INTEGER  TYPE, I, J 

REAL*8  TO(0:M,0:N) ,TN(0:M,0:N) , ALPHA(NT) , RHOC (NT) ,DR,DT,C1 
REAL*8  R0,R1,C2(NT),C,QM(0:M,0:N) 

COMMON  TO,TN,QM 
R1=R0+DR 

Cl=2.0*ALPHA(TYPE)*DT/( (Rl**2-R0**2 )*DR) 

C=1 . 0-Cl*R0-Cl*Rl-2 . 0*C2 ( TYPE ) 

IF  (C.LT.0.0)  THEN 

PRINT*, ’INTR  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN( I , J )=C1*(RO*TO( I , J-l )+Rl*TO( I , J+l ) )+C2(TYPE)* 

.  (TO( 1-1 , J )+TO( 1+1 , J) )+C*TO( I , J )+QM( I , J )*DT /RHOC (TYPE) 

RETURN 
END 


SUBROUTINE  TO  CALCULATE  TPRIME  (STORED  IN  T)  AT  DT  FOR  A  NODE 
AT  A  BOUNDARY  LAYER  WITH  TWO  MATERIAL  TYPES. 

SUBROUTINE  R2MATL ( K , RHOC , DX , DR , DT , TYP1 , TYP2 , I , J , RO ) 

INCLUDE  ’ head . h ’ 

INTEGER  I ,J,TYP1,TYP2 

REAL*8  TO(0:M,0: N) ,TN(0:M,0:N) ,K(NT) , RHOC (NT) ,DX,DR,DT 
REAL*8  C , Cl , C2 , C3 , C4 , CXX , RO , R1 , R2 , QM ( 0 : M , 0 : N ) 

COMMON  TO,TN,QM 

Rl=R0+DR/2.0 

R2=R0+DR 

CXX=RHOC(TYP1)*(R1**2-RO**2)+RHOC(TYP2)*(R2**2-R1**2) 

Cl =2 . 0*K ( TYP1 ) *DT/ ( CXX*DR ) 

C2=2 . 0*K(  TYP2 )  *DT/  ( CXX*DR ) 

C3=K(TYP1 )*(R1**2-R0**2)*DT/(CXX*DX**2 ) 

C4=K(TYP2 )*(R2**2-R1**2 )*DT/(CXX*DX**2 ) 

C=1 . 0-Cl*R0-C2*R2-2 . 0*C3-2 . 0*C4 
IF  (C.LT.0.0)  THEN 

PRINT*,  ’R2MATL  UNSTABLE  CONDITION;  C=\C 
ENDIF 

TN( I , J )=C1*RO*TO( I , J-l )+C2*R2*TO( I , J+l )+C3*(TO( 1-1 , J ) 

+TO( 1+1 , J ) }+C4*(TO(I-l , J )+TO(I+l , J) )+QM( I , J )*(R2**2 
-R0**2 ) *DT/CXX+C*TO( I , J ) 

RETURN 

END 

SUBROUTINE  TO  ACCOUNT  FOR  INSULATIVE  SIDE  BOUNDARY  CONDITIONS  FOR  A 
SINGE  ELEMENT  NODE. 

SUBROUTINE  SIDE1(ALPHA, RHOC, TYPE, DR, DT, I , J,C2,R0, IFIN) 

INCLUDE  ’head. h’ 

INTEGER  I, J, TYPE, IFIN 

REAL*8  TO(0:M,0:N),TN(0:M,0:N) , ALPHA (NT) , RHOC (TYPE) ,DR,DT 
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REAL*8  C,C1,C2(NT),R0,R1,QM(0:M,0:N) 

COMMON  TO.TN.QM 
R1=R0+DR 

Cl=2 . 0* ALPHA ( TYPE ) *DT/ ( ( Rl**2 -R0**2 ) *DR ) 

C=1 . 0-Cl*Rl-Cl*R0-2 . 0*C2 ( TYPE ) 

IF  (C.LT.0.0)  THEN 

PRINT*,  ’SIDE1  UNSTABLE  CONDITION;  C=\C 
ENDIF 

IF  (I.EQ.O)  K=I+1 
IF  (I.EQ.IFIN)  K=I-1 

TN(I , J)=C1*R1*T0( I , J+1)+C1*RO*TO(I , J-l )+C2 (TYPE) *2.0* 

.  TO(K, J)+QM( I , J )*DT/RHOC(TYPE)+C*TO( I , J ) 

RETURN 

END 

SUBROUTINE  TO  ACCOUNT  FOR  THE  INSULATIVE  SIDE  BOUNDARY  CONDITIONS 
FOR  A  NODE  CONTANING  TWO  DIFFERENT  MATERIALS 

SUBROUTINE  SIDE2(K,RH0C,DX,DR,DT,TYP1,TYP2,I ,J,RO,IFIN) 

INCLUDE  ’head.h’ 

INTEGER  TYP1 , TYP2 , 1 , J , KK , I  FI N 

REAL*8  TO(0:M,0:N),TN(0:M,0:N) ,RHOC(NT) ,K(NT) ,DX ,DR,DT 
REAL*8  CXX , C , Cl , C2 , C3 , C4 , RO , R1 , R2 , QM ( 0 : M , 0 : N ) 

COMMON  TO,TN,QM 

Rl=R0+DR/2.0 

R2=R0+DR 

CXX=RHOC(TYP1)*(R1**2-RO**2)+RHOC(TYP2)*(R2**2-R1**2) 

C1=2.0*K(TYP1)*DT/(CXX*DR) 

C2=2.0*K( TYP2 ) *DT/ ( CXX*DR ) 

C3=K(TYP1)*(R1**2-R0**2)*DT/(CXX*DX**2) 

C4=K(TYP2)*(R2**2-R1**2)*DT/(CXX*DX**2) 

IF  (I.EQ.O)  KK=I+1 

IF  (I.EQ.IFIN)  KK=I-1 

C=1 . 0-Cl*R0-C2*R2-2 . 0*C3-2 . 0*C4 

IF  (C.LT.0.0)  THEN 

PRINT*, ’SIDE2  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I ,J)=C1*RO*TO(I ,J-l)+C2*R2*TO(I , J+l )+2.0*C3* 

TO(KK,J)+2.0*C4*TO(KK,J)+QM(I,J)*(R2**2-RO**2)*DT 
.  /CXX+C*TO( I , J) 

RETURN 

END 

SUBROUTINE  TO  CALCULATE  THE  TEMPERATURE  EFFECTS  AT  THE  INSULATED 
CORNERS  REPRESENTING  THE  SEMI-INFINITE  BOUNDARY  CONDITIONS  OF 
THE  MUSCLE  TISSUE. 

SUBROUTI NE  I NSCNR (ALPHA , RHOC , DR , DT , I , J , TYPE , C2 , RO , I FI N ) 

INCLUDE  * head.h’ 

INTEGER  I ,J, TYPE, I  FIN 

REAL*8  TO(0:M,0:N),TN(0:M,0:N) , ALPHA ( NT ), RHOC (NT) 

REAL*8  C,C1 ,C2(NT) ,R0,R1 ,DR,DT,QM(0:M,0:N) 

COMMON  TO.TN.QM 
Rl=R0+DR/2.0 
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Cls2.0*ALPHA(TYPE)*DT/((Rl**2-R0**2)*DR) 

C=1.0-C1*R0-2.0*C2(TYPE) 

IF  (C.LT.O.O)  THEN 

PRINT*, ’INSCNR  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

IF  (I.EQ.O)  K=I+1 
IF  (I.EQ.IFIN)  K=I-1 

TN(I,J)=C1*RO*TO(I,J-1)+2.0*C2(TYPE)*TO(K,J)+QM(I,J)* 

.  DT/RHOC ( TYPE ) +C*TO ( I , J ) 

RETURN 

END 

SUBROUTINE  TO  ACCOUNT  FOR  THE  INSULATIVE  MUSCLE  LAYER 

SUBROUTINE  I NSBND ( ALPHA , RHOC , DR , DT , TYPE , I , J , C2 , RO ) 

INCLUDE  ’head.h’ 

INTEGER  TYPE, I, J 

REAL*8  TO(0:M,0:N),TN(0:M,0:N), ALPHA ( NT ) , RHOC ( NT ) , DR , DT 
REAL*8  C,C1,C2(NT) ,R0,R1 ,QM(0:M,0:N) 

COMMON  TO,TN,QM 
Rl=R0+DR/2.0 

Cl =2 . 0* ALPHA ( TYPE ) *DT/ ( ( R1 **2-R0**2 ) *DR ) 

C=1 . 0-Cl*R0-2 . 0*C2 (TYPE ) 

IF  (C.LT.O.O)  THEN 

PRINT*,’ I  NSBND  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I ,J)=C1*RO*TO(I ,J-1)+C2(TYPE)*(T0(I-1,J)+T0(I+1,J) )+ 

C*TO { I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 

RETURN 

END 

SUBROUTINE  TO  ACCOUNT  FOR  THE  INSULATIVE  CORNER  BOUNDARY  CONDITIONS 
WHERE  CONVECTION  OCCURS  AT  THE  BLOOD/PLAQUE  INTERFACE 

SUBROUTINE  CONVCNR ( H , ALPHA , RHOC , DR , DT , TYPE , 

I , J ,TB,C2 ,R0 , IFIN) 

INCLUDE  ’head.h’ 

INTEGER  TYPE, I, J,K, IFIN 

REAL*8  TO(0:M,0:N) ,TN(0:M,0:N) ,ALPHA(NT) , RHOC (NT) ,DR,DT 
REAL*8  H(0:M) ,C,C1 ,C2(NT) ,C3,R1 ,R0,TB(0:M+2) ,QM(0:M,0:N) 

COMMON  TO,TN,QM 
Rl=R0+DR/2.0 
IF  (I.EQ.O)  K=I+1 
IF  (I.EQ.IFIN)  K=I-1 

Cl=2. 0*ALPHA(TYPE)*DT/( (R1**2-R0**2)*DR) 

C3=2.0*H(I )*RO*DT/(RHOC(TYPE)*(R1**2-RO**2) ) 

C=1.0-C1*R1-2.0*C2(TYPE)-C3 

IF  (C.LT.O.O)  THEN 

PRINT*, ’CONVCNR  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I , J)=C1*R1*T0( I , J+l )+2.0*C2(TYPE)*TO(K, J)+C*TO( I , J ) 

.  +C3*TB( 1+1 )+{QM( I , J )/RHOC(TYPE) )*DT 

RETURN 
END 
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SUBROUTINE  TO  CALCULATE  TPBIME  AT  DT  FOR  A  RODE  LYING  ON  A 
CONVECTIVE  BOUNDARY  LAYER  (PLAQUE/BLOOD) 

SUBROUT I NE  CONVBND ( H , ALPHA , RHOC , DR , DT , TYPE , I , J , TB , C2 , RO ) 

INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) , ALPHA ( NT ), RHOC (NT) 

.  , QM ( 0 : M , 0 : N ) , C2 ( NT ) , C , DR , DT , C 1 , C3 , RO , R 1 , TB ( 0 : M+ 2 ) 

COMMON  TO , TN , QM 
Rl=R0+DR/2.0 

Cl=2 . 0*ALPHA( TYPE ) *DT/( ( Rl**2-R0**2 )*DR ) 
C3=2.0*H(I)*R0*DT/(RHOC(TYPE)*(Rl**2-R0**2)) 

C=1 . 0-Cl*Rl-2 . 0*C2 ( TYPE ) -C3 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONVBND  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I , J)=C1*R1*T0( I , J+l)+C2(TYPE)*(TO( 1-1 , J )+TO(I+l, J) ) 

+C3*TB ( I + 1 ) +C*TO ( I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 

RETURN 

END 

SUBROUTINES  FOR  THE  VARIOUS  CONVECTIVE  BOUNDARY  LAYER  CONDITIONS  IN 
THE  PLAQUE  LAYER  OF  VARYING  RADII. 

SUBROUTI NE  CONV3A ( H , ALPHA , RHOC , DR , DX , DT , TYPE , I , J , TB , CC , RO ) 
INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) , ALPHA ( NT ), RHOC (NT) , 

.  QM(0:M,0:N) ,C2,C,DR,DT,C1 ,C3,R0,R1 ,TB(0:M+2) ,CC(NT) ,DX 

COMMON  TO,TN,QM 
R1=R0+DR 
C1=2.0*CC(TYPE) 

C2=2. 0*ALPHA ( TYPE ) *DT/ { DR* ( Rl**  2 -R0**2 ) ) 

C3=2 . 0*H ( I ) *DT/ ( RHOC ( TYPE ) *DX ) 

C=1 . 0-Cl-C2*Rl-C2*R0-C3 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONV3A  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN( I , J )=Cl*TO( 1-1 , J )+C2*Rl*TO( I , J+l )+C2*RO*TO( I , J-l) 

+C3*TB( I +2 ) +C*TO ( I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 

RETURN 

END 


SUBROUTINE  CONV3B(H , ALPHA , RHOC , DR ,DX , DT , TYPE , I , J , TB, CC , RO ) 
INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL*8  TO(0:M,0:N) ,TN(0:M,0:N),H(0:M) , ALPHA ( NT ), RHOC (NT) , 

.  QM(0:M,0:N) ,C2,C,DR,DT,C1 ,C3,R0,R1 ,TB(0:M+2) ,CC(NT) ,DX 

COMMON  TO,TN,QM 
R1=R0+DR 
ri=2.0*CC(TYPE) 

C2=2 . 0*ALPHA ( TYPE ) *DT/ ( DR* ( Rl **2-R0**2 ) ) 
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C3=2 . 0*H ( I ) *DT/ ( BHOC ( TYPE )  *DX ) 

C=1 . 0-Cl-C2*Rl-C2*R0-C3 
IF  (C.LT.0.0)  THEN 

PRINT*,’ CONV3B  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I,J)=C1*TO(I+1,J)+C2*R1*TO(I,J+1)+C2*RO*TO(I,J-1) 
+C3*TB( I ) +C*TO ( I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 

RETURN 

END 


SUBROUTINE  CONV4A(H,ALPHA,RHOC,DR,DX,DT,TYPE,I ,J ,TB,CC,R0) 
INCLUDE  ’head.h’ 

INTEGER  I ,J,TYPE 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M), ALPHA ( NT ) , RHOC ( NT ) , DX 
.  ,QM(0:M,0: N) ,C2,C,DR,DT,C1 ,C3,C4,R0,Rl,TB(0:M+2) ,CC(NT) ,CR 

COMMON  TO,TN,QM 
Rl=R0+DR/2.0 
CR=R1**2-R0**2 
C1=2.0*CC(TYPE) 

C2=2. 0*ALPHA ( TYPE ) *DT*R1/ (DR*CR) 

C3=2 . 0*H ( I ) *DT/ ( RHOC ( TYPE ) *DX ) 

C4=2 . 0*H ( I ) *R0*DT/ ( RHOC ( TYPE ) *CR ) 

C=1 . 0-C1-C2-C3-C4 
IF  (C.LT.0.0)  THEN 

PRINT*, 'CONV4A  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN( I , J )=Cl*TO( 1+1 , J )+C2*TO( I , J+l )+C3*TB( I )+C4*TB( 1+1 )+C* 

.  TO( I , J)+QM( I , J)*DT/RHOC(TYPE) 

RETURN 

END 


SUBROUTINE  OONV4B ( H , ALPHA , RHOC , DR , DX , DT , TYPE . I , J , TB , CC , RO ) 
INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) , ALPHA (NT) , RHOC (NT) ,DX 
.  ,QM(0:M,0:N) ,C2,C,DR,DT,C1 ,C3,C4,R0,R1 ,TB(0:M+2) ,CC(NT) ,CR 

COMMON  TO,TN,QM 
Rl=R0+DR/2.0 
CR=R1**2-R0**2 
C1=2.0*CC(TYPE) 

C2=2 . 0* ALPHA ( TYPE ) *DT*R1/ ( DR*CR ) 

C3=2 . 0*H ( I ) *DT/ ( RHOC ( TYPE ) *DX ) 

C4 = 2 . 0*H ( I ) *R0*DT/ ( RHOC ( TYPE ) *CR ) 

C= 1 . 0-C1-C2-C3-C4 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONV4B  UNSTABLE  CONDITION;  C*’,C 
ENDIF 

TN(I , J)=C1*T0(I-1 , J)+C2*TO( I , J+l )+C3*TB( I+2)+C4*TB(I tl )+C* 

.  TO( I , J )+QM( I , J )*DT/RHOC(TYPE) 

RETURN 

END 
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SUBROUTINE  CONV5A(H,ALPHA,RHOC,DR,DX,DT,TYPE,I , J,TB,R0) 

INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) , ALPHA (NT) ,RHOC(NT) , 

.  C2 , DR , DT , Cl , C3 , C4 , RO , R1 , R2 , TB ( 0 : M+2 ) , C5 , C6 , VOL , PI , CR2 ,  DX ,  C 

.  ,QM(0:M,0:N),CR1 

COMMON  TO,TN,QM 
PI=3. 141592654 
Rl=R0+DR/2 . 0 
R2=RO+DR 
CR1=R2**2-R1**2 
CR2=R2**2-R 0**2 
VOL=PI * ( CR2*DX-CR1 *DX/2 . 0 ) 

C1=ALPHA ( TYPE ) *PI *CR2*DT/ ( VOL*DX ) 

C2=ALPHA ( TYPE )  *PI *CR1*DT/ ( VOL*DX ) 

C3=ALPHA ( TYPE )  *PI *R0*DX*DT/ ( VOL*DR ) 

C4=ALPHA ( TYPE ) *2 . 0*PI *R2*DX*DT/( VOL*DR ) 

C5=H( I )*PI *R1*DX*DT/ ( VOL*RHOC ( TYPE ) ) 

C6=H ( I )*PI *CR1*DT/ ( VOL*RHOC ( TYPE ) ) 

C=1 . 0-C1-C2-C3-C4-C5-C6 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONV5A  UNSTABLE  CONDITION;  C=*,C 
ENDIF 

TN( I , J)=Cl*TO( 1-1, J )+C2*TO( 1+1 , J )+C3*TO( I , J-l )+C4*TO( I , J+l ) 

.  +C5*TB ( I +1 ) +C6*TB( I +2 ) +C*TO( I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 

RETURN 
END 

SUBROUTI NE  CONV5B(H , ALPHA , RHOC , DR , DX , DT , TYPE , I , J , TB , RO ) 

INCLUDE  ’head.h’ 

INTEGER  I, J, TYPE 

REAL* 8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) , ALPHA ( NT ), RHOC (NT) , 

.  C2 , DR , DT , C 1 , C3 , C4 , RO , R1 , R2 , TB ( 0 : M+2 ) , C5 , C6 , VOL , PI , CR2 , DX , C 

,QM(0:M,0:N),CR1 
COMMON  TO,TN,QM 
PI=3. 141592654 
Rl=R0+DR/2.0 
R2=RO+DR 
CR1=R2**2-R1**2 
CR2=R2**2-R0**2 

VOL=PI*(CR2*DX-(Rl**2-RO**2)*DX/2. 0) 

C1=ALPHA ( TYPE ) *PI *CR2*DT/ ( VOL*DX ) 

C2=ALPHA ( TYPE ) *PI *CR1*DT/ ( VOL*DX ) 

C3=ALPHA ( TYPE ) *PI *R0*DX*DT/( VOL*DR ) 

C4= ALPH A ( TYPE ) *  2 . 0*PI *R2  *DX*DT/ ( VOL*DR ) 

C5=H(I )*PI*Rl*DX*DT/(VOL*RHOC(TYPE) ) 

C6=H(I )*PI*(R1**2-RO**2)*DT/(VOL*RHOC(TYPE) ) 

C=1 . 0-C1-C2-C3-C4-C5-C6 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONV5B  UNSTABLE  CONDITION;  C=\C 
ENDIF 

TN ( I , J)=C1*T0(I+1, J )+C2*TO( 1-1 , J)+C3*TO( I , J-l )+C4*TO( I , J+l ) 
+C5*TB ( I + 1 ) +C6*TB ( I ) +C*TO ( I , J ) +QM ( I , J ) *DT/RHOC ( TYPE ) 
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RETURN 

END 

SUBROUTINE  CONV6A(H,K,RHOC,DR,DX,DT,TYP1,TYP2,I,J,TB,RO) 
INCLUDE  ’head.h* 

INTEGER  TYP1,TYP2,I ,J 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M) ,K(NT) ,RHOC(NT) , 

.  DR,DX,DT,RO,R1,R2,C,VOL,PI,CR1,CR2,C1,C2,C3,C4,C5,C6,C7, 

.  QM(0:M,0:N} ,TB(0:M+2) 

COMMON  TO,TN,QM 

PI=3. 141592654 

Rl=R0+DR/2.0 

R2=R0+DR 

CR1=R2**2-R1**2 

CR2=R1**2-R0**2 

VOL=PI * ( ( R2**2-R0**2 ) *DX-CR2*DX/2 . 0 ) 

Cl = ( PI * ( RHOC ( TYP1 ) *CR2  *DX /2 . 0+RHOC ( TYP2 ) *CR1 *DX ) ) /DT 
C2=K ( TYP1 ) *PI *CR2/ ( DX  *C 1 ) 

C3=K(TYP2)*PI*CR1/(DX*C1) 

C4=K ( TYP1 ) *PI *R0*DX/ ( DR*C1 ) 

C5=2 . 0*K ( TYP2 ) *PI *R2*DX/ ( DR*C1 ) 

C6=H(I)*PI*R1*DX/C1 
C7=H( I )*PI*CR2/C1 
C=1 . 0-C2-2 . 0*C3-C4-C5-C6-C7 
IF  (C.LT.0.0)  THEN 

PRINT*, ’CONV6A  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(I , J)=C2*TO(I-l , J)+C3*(TO(I-l , J)+TO( 1+1 , J) )+C4*TO( I , J-l) 
.+C5*TO( I , J+l)+C6*TB(I+l)+C7*TB(I+2)+C*TO(I , J)+QM(I , J)*V0L/C1 
RETURN 
END 


SUBROUTI NE  C0NV6B (H,K, RHOC , DR , DX , DT , TYP1 , TYP2 , I , J , TB , RO ) 
INCLUDE  ’head.h* 

INTEGER  TYP1,TYP2,I ,J 

REAL*8  TO(0:M,0:N),TN(0:M,0:N),H(0:M),K(NT) ,RHOC(NT) , 

.  DR , DX , DT , RO , R1 , R2 , C , VOL , PI , CR1 , CR2 , Cl , C2 , C3 , C4 , C5 , C6 , C7 , 

.  QM(0:M,0:N) ,TB(0:M+2) 

COMMON  TO,TN,QM 
PI=3. 141592654 
Rl=R0+DR/2 . 0 
R2=R0+DR 
CR1=R2**2-R1**2 
CR2=R1**2-R0**2 

VOL=PI * ( ( R2**2-RO**2 ) *DX-CR2*DX/2 . 0 ) 

Cl= ( PI * ( RHOC( TYP1 ) *CR2*DX/2 . O+RHOC ( TYP2 ) *CR1*DX ) ) /DT 
C2=K(TYP1 )*PI*CR2/(DX*C1 ) 

C3=K ( TYP2 ) *PI *CR1/ ( DX*C1 ) 

C4=K(TYP1 )*PI*R0*DX/(DR*C1 ) 

C5=2 . 0*K ( TYP2 ) *PI *R2*DX/( DR*C1 ) 

C6=H(I)*PI*R1*DX/C1 

C7=H(I)*PI*CR2*DX/C1 

C=1 . 0-C2-2 . 0*C3-C4-C5-C6-C7 
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IF  (C.LT.0.0)  THEN 

PRINT*, '(X)NV6B  UNSTABLE  CONDITION;  C=’,C 
ENDIF 

TN(  I , J)=C2*T0(I+1 , J)+C3*(TO(I+l , J)+TO( 1-1 , J) )+C4*TO(I , J-l ) 
.  +C5*TO (I , J+l) +C6*TB( I + 1 ) +C7*TB( I ) +C*TO ( I , J ) +QM ( I , J ) *VOL/Cl 
RETURN 
END 


DOUBLE  PRECISION  FUNCTION  COEFF( ALPHA, DT,DX, TYPE) 
include  ’head.h’ 

INTEGER  TYPE 

REAL*8  ALPHA ( NT ),DT,DX 

COEFF=ALPHA ( TYPE ) *DT/DX**2 

RETURN 

END 


CUBIC  SPLINE  INTERPOLATION  SUBROUTINE  FOR  CALCULATING  THE  AVERAGE 
CONVECTIVE  HEAT  TRANSFER  COEFFICIENTS  AND  ENERGY  DEPOSITION  FROM  THE 
MICROWAVE  APPLICATOR  AT  THE  NODE  INTERFACES. 

SUBROUTINE  CUBSPL ( X , FX , F , XBAR , NN , NX ) 

INCLUDE  ’head.h’ 

INTEGER  NN,NX 
PARAMETER  (N0=200) 

REAL*8  X(0:NN) ,F(0:NN) ,XBAR(0:NX) ,FX(0:NX) 

REAL*8  ALPHA (NO) ,B(0:NO) ,C(0:NO) ,D(0:NO) ,H(0:NO) 

,L(0:NO) ,MU(0:NO) ,Z(0:NO) ,X1 
DO  10  I =0 , NN-1 
H ( I )=X(I+1)-X(I ) 

10  CONTINUE 

DO  20  1=1, NN-1 

ALPHA(I)=3.0*(F(I+1)*H(I-1)-F(I)*(X(I+1)-X(I-1))+F(I-1)*H(I)) 
/(H( 1-1 )*H( I ) ) 

20  CONTINUE 
L(0)=1.0 
MU ( 0 ) =0 . 0 
Z(0)=0.0 
DO  30  1=1, NN-1 

L(I)=2.0*(X(I+1)-X(I-1) )-H(I-l)*MU(I-l) 

MU(I)=H(I)/L(I) 

Z(I )=(ALPHA(I )-H(I-l)*Z(I-l) )/L(I ) 

30  CONTINUE 
L(NN)=1,0 
Z(NN)=0.0 
C(NN)=0.0 
DO  40  I=NN-1 ,0,-1 

C(I)=Z(I)-MU(I)*C(I+1) 

B(I)=(F(I+1)-F(I))/H(I)-H(I)*(C(I+1 )+2.0*C(  I ) )/3.0 

D(I)=(C(I+1)-C(I) )/(3.0*H(  I ) ) 

40  CONTINUE 
DO  50  1=0, NX 
J=0 
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IF  (XRAR(I).LT.X(O))  THEN 
FX ( I )=F(0) 

GO  TO  50 
ENDIF 

IF  (XBAR(I).GT.X(NN))  THEN 
FX(I)=F(NN) 

GO  TO  50 
ENDIF 

DO  55  K=1 ,NN 

IF  ( (XBAR( I ) . LT. X(K) ) . AND* (XBAR( I ) .GE. X(K-1 ) ) )  J=K-1 
55  CONTINUE 

IF  (F(J).EQ.F(J+1))  THEN 
FX( I )=F( J ) 

GO  TO  50 
ENDIF 

X1=XBAR( I )-X( J ) 

FX(I)=F(J)+B(J)*X1+C(J)*(X1**2)+D(J)*(X1**3) 

50  CONTINUE 
RETURN 
END 


C 

C  SUBROUTINE  TO  ESTIMATE  THE  CONVECTIVE  HEAT  TRANSFER  COEFFICIENTS 
C  OF  THE  BLOOD  VESSEL  WALL  AND  MICROWAVE  APPLICATOR  USING  A  CUBIC 
C  SPLINE  APPROXIMATION  FOR  THE  ANALYTICAL  SOLUTION  OF  A  CONCENTRIC 
C  ANNULAR  CIRCULAR  PIPE. 

C 

SUBROUTI NE  HTC ( RPLX , RPRX , I FI N , HP, HM , K ) 

INCLUDE  ’head.h* 

INTEGER  I FIN, I 
PARAMETER  ( N0=200 , TYPE=5 ) 

REAL*8  RSTR1(0:6) ,RSTR2{0: 7 ) ,NU1(0:6) ,NU2(0: 7) ,RPRX(0:M) , 

.  NU(0:NO),RBAR(0:NO),K(NT),HP(0:M),HM(0:M),DH,RPLX(0:M) 

DATA  (RSTR1( I), 1=0, 6)/. 05,. 10, .20,. 40,. 60,. 80, 1.00/ 

DATA  (NU1( I), 1=0, 6)/17. 81, 11. 91, 8. 499, 6. 583, 5. 912, 5. 58, 5. 385/ 
DATA  (RSTR2(I ) , 1=0, 7)/0.0, .05 , . 10, . 20, . 40, . 60, .80,1.00/ 

DATA  (NU2 (I), 1=0, 7) /4. 364, 4. 792, 4. 834, 4. 883 ,4.979,5.099, 
5.24,5.385/ 

DO  10  I=0,IFIN 

RBAR(I )=RPRX(I )/RPLX(I ) 

10  CONTINUE 

CALL  CUBSPL ( RSTR1 , NU , NU 1 , RBAR , 6 , I FI N ) 

DO  20  1=0, I FIN 
IF  (RPRX(I).EQ.O.O)  THEN 
HM(I )=0.0 
ELSE 

DH=2 . 0* ( RPLX ( I ) -RPRX ( I ) } 

HM( I )=NU( I )*K(TYPE)/DH 
ENDIF 

20  CONTINUE 

CALL  CUBSPL ( RSTR2 , NU , NU2 , RBAR , 7 , I FI N ) 

DO  30  1=0, IFIN 
IF  (RPRX(I).EQ.O.O)  THEN 
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DH=2.0*RPLXd) 

NU=3.66 

HP(I)=NU*K(TYPE)/DH 

ELSE 

DH=2 . 0* { RPLX ( I ) -RPRX (I ) ) 

HP( I )=NU(I )*K(TYPE)/DH 
C  PRINT*, ’H,NU  1  ’.HPdJ.NUd) 

END  IF 

30  CONTINUE 
RETURN 
END 

SUBROUTINE  TO  CALCULATE  THE  HEAT  DISTRIBUTION  IN  THE  BLOOD  LAYER 
DUE  TO  ONE-DIMENSIONAL  CONVECTIVE  HEAT  TRANSFER 

SUBROUTI NE  BFLOW ( TB1 , QB , RHOC , TYPE , HPL , HPR , RPL , RPR , 
DX,DT,TP,V,IFIN, JRPR , JRPL , QAVE ) 

INCLUDE  ’head.h’ 

INTEGER  TYPE, IFIN,L, JRPR(0:M) , JRPL(0:M) 

REAL*8  TN(0:M,0:N) ,TBl(0:M+2 ) ,RPL(0:M) ,RPR(0:M) ,QB(0:M) , 

.  HPL(0:M) ,HPR(0:M) ,V(0:M) ,DX,DT,C,C1,C2,C3,C4,TB2( 0:452) 

.  ,RHOC(NT),TP(0:M),QAVE,QM(0:M,0:N),TO(0:M,0:N) 

COMMON  TO,TN,QM 
DO  10  1=0, IFIN 

Cl=2. 0*DT/ ( RHOC ( TYPE ) * ( RPL (I ) **2-RPR (I ) **2 ) ) 

C2=HPR( I )*RPR( I )*C1 
C3=HPL( I )*RPL( I )*C1 

V(I )=QAVE*5 . 3051647E-9/(RPL( I )**2-RPR(I )**2 ) 

C4=V ( I )*DT/DX 
C=1 . 0-2 . 0*C4-C2-C3 
IF  (C.LT.0.0)  THEN 

PRINT*, ’BFLOW  UNSTABLE  CONDITION,;  C=*,C 
ENDIF 

TB2 (I + 1 ) =C4* ( TB1 (I ) +TB1 (I + 2 ) ) +C2  *TP (I ) +C3*TN (I , JRPL (I ) ) +C* 

.  TB1 ( I + 1 ) +QB ( I ) *DT/RHOC ( TYPE ) 

10  CONTINUE 

TB1 ( I  FI N+2 ) =TB2 ( I  FI N+ 1 ) 

DO  20  1=1 , IFIN+1 
TB1 ( I )=TB2( I ) 

20  CONTINUE 

DO  40  1=0, IFIN 

DO  30  J=JRPR( I ) , JRPL( I ) — 1 
TN( I , J )=TB1(I+1 ) 

30  CONTINUE 
40  CONTINUE 
RETURN 
END 

SUBROUTINE  TO  RECALCULATE  THE  POWER  GENERATED  BY  THE  MICROWAVE 
APPLICATOR  AT  THE  INTERFACE  OF  THE  NODES  USING  THE  CUBIC  SPLINE 
SUBROUTINE. 

SUBROUT I NE  I NTRP1 ( I FLAG , AOUT , DR , DX ) 

INCLUDE  ’head.h’ 


I NTEGER  I , J , I START ,  I  FI N , J START , JFI N , MODEL (M , N ) , NN1 , NN2 
INTEGER  NX 1, NX 2 

REAL*8  A0UT(0:M,0:N) ,X(0:200) ,F( 0:200) , FX(0: 200) ,DR,DX 
REAL*8  XBAR(0: 200) ,RBAR(0: 200) ,TEMP(M,N) 

REAL*8  POWER(M,N) ,QM(0:M,0:N) ,TO{0:M,0:N) ,TN(0:M,0:N) 

COMMON  TO.TN.QM 

COMMON  /DAT/POWER , TEMP , MODEL , I START , I  FI N , JSTART , JFI N 

NN 1=1 FIN- I START 

NN2=JFIN-JSTART 

NX1=IFIN 

NX2=JFIN 

DO  10  J=JSTART, JFIN 
DO  20  I =1 START, I  FIN 
IF  (IFLAG.EQ.l)  THEN 
F( 1-1 )=POWEB(I , J ) 

ELSE 

F(I-1)=TEMP(I,J) 

ENDIF 

X(I-l)=I*DX-DX/2.0 
XBAR( 1-1 )=( 1-1 )*DX 
20  CONTINUE 

XBAR( IFIN)=XBAR( IFIN-1 )+DX 
CALL  CUBSPL(X,FX,F,XBAR,NN1,NX1) 

DO  30  I=ISTART-1 , IFIN 
AOUT( I , J )=FX( I ) 

30  CONTINUE 
10  CONTINUE 

DO  40  I =1 START-1, 1 FIN 
DO  50  J=JST ART, JFIN 
F( J-l)=AOUT(I , J) 

X( J-l)=J*DR-DR/2.0 
RBAR( J-l )=( J-1)*DR 
50  CONTINUE 

RBAR( JFIN)=RBAR(JFIN-1 )+DR 
CALL  CUBSPL ( X , FX , F , RBAR , NN2 , NX2 ) 

DO  60  J=JSTART-1 , JFIN 
AOUT(I ,J)=FX(J) 

60  CONTINUE 

40  CONTINUE 
RETURN 
END 
C 

C  ALTERNATIVE  INTERPOLATION  SUBROUTINE  FOR  COURSE  MESH  ANGIOPLASTY 
C  MODELS  WHERE  LARGE  DIFFERENCES  IN  SAR  VALUES  OCCUR  IN  THE  EE  DATA 
C  GENERATED.  SUCH  LARGE  DIFFERENCES  DO  NOT  PERMIT  ACCURATE 
C  APPLICATION  OF  THE  CUBIC  SPLINE  SUBROUTINE 
C 

SUBROUTI NE  I NTRP2 ( I  FLAG , AOUT , DR , DX ) 

INCLUDE  ’head.h’ 

INTEGER  I ,J , ISTART, IFIN, JFIN, MODEL(M,N) 

REAL*8  AOUT(0:M,0:N) ,DR,DX,TEMP(M,N) , POWER (M,N) ,QM(0:M,0:N) 
REAL*8  TO(0:M,0:N),TN(0:M,0:N),'AIN(M,N) 

COMMON  TO.TN.QM 

COMMON  /DAT/POWER , TEMP , MODEL , I START , I FI N , JSTART .JFIN 
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DO  10  J=JSTART, JFIN 
DO  20  I=ISTART , IFIN 
IF  (IFLAG.EQ.l)  THEN 
AIN(I,J)=POWER(I,J) 

ELSE 

AIN(I,J)=TEMP(I,J) 

END  IF 

20  CONTINUE 

10  CONTINUE 

DO  30  I=ISTART, IFIN-1 
DO  40  J=JSTART, JFIN-1 

AOUT( I ,J)=(AIN(I,J )+AIN(I +1 , J )+AIN( I ,  J+l )+AIN( 1+1 , J+l ) )/4 . 0 
40  CONTINUE 

30  CONTINUE 

DO  50  1=0, IFIN, IFIN 
IF  (I.EQ.O)  L=1 
IF  (I.EQ.IFIN)  L=IFIN 
DO  60  J=JSTART, JFIN 
AOUT(I ,  J)  =  (AIN(L,  J)+.\IN(L,  J+l )  )/2.0 
60  CONTINUE 
50  CONTINUE 

DO  70  J=0, JFIN, JFIN 
IF  (J.EQ.O)  L=1 
IF  (J.EQ.JFIN)  L=JFIN 
DO  80  I=IST ART, IFIN 

AOUT( I,J)  =  (AIN(I , L)+AIN( 1+1 , L) )/2. 0 
80  CONTINUE 

70  CONTINUE 

AOUT ( 0 , 0 ) = ( AOUT (0,1) +AOUT ( 1 , 0 ) ) /2 . 0 

AOUT( IFIN, JFIN)=( AOUT( IFIN-1, JFIN)+AOUT{ IFIN, JFIN-1 ))/2.0 
AOUT ( 0 , J  FI N )  =  ( AOUT ( 0 , J  FI N- 1 ) + AOUT (l,JFIN))/2.0 
AOUT( IFIN, 0)=( AOUT( IFIN-1, 0)+AOUT( IFIN, l))/2.0 
999  RETURN 
END 

SUBROUTINES  PDATAl  AND  PDATA2  FORMAT  DATA  OUTPUT 

SUBROUTINE  PDATA 1 ( I FN , OUTP1 , T I TLE , T I ME ) 

INCLUDE  ’head.h’ 

PARAMETER  (LINES=17,COLM=12) 

INTEGER  I , J , I START, IFIN, JSTART, JFIN, MODEL (M,N) , K, 1 1 , 12 , IFN 
REAL*8  POWER (M,N) ,QM(0:M,0:N),TO(0:M,0:N),TN(0:M,0:N) 

, TEMPI M,N) ,OUTP1(0:M,0:N) ,TIME 
CHARACTER* 40  TITLE 

COMMON/DAT/POWER , TEMP , MODEL , I START , I FI N , JSTART , JFIN 
COMMON  TO,TN,QM 
WRITE  (IFN, 300)  TIME, TITLE 
300  FORMAT  (IX, *TIME=’ , E10.3,2X,40A) 

DO  100  K=l, LINES 
Il=COLM*(K-l) 

I2=I1+C0LM-1 
IF  (K.EQ. LINES)  I2=IFIN 
WRITE  (IFN, 400)  (1,1=11,12) 

FORMAT  (IX, 12  110) 


400 


DO  200  J=JSTART-1, JFIN 

WRITE  (IFN.500)  J, (0UTP1( I , J) , 1=1 1, 12) 

500  F0RMAT(1X,I3,1X,12  E10.3) 

200  CONTINUE 
100  CONTINUE 
RETURN 
END 

SUBROUTINE  PDATA2 ( I FN , OUTP2 , TI TLE ) 

INCLUDE  'head.h* 

PARAMETER  (LINES=17,COLM=12) 

INTEGER  I ,J, I START, IFIN, JSTART, JFIN, MODEL (M,N) , K, 1 1 , 12 , 1 FN 
REAL*8  POWER(M,N) ,QM(0:M,0:N) ,TO(0:M,0:N),TN(0:M,Q:N) 
,TEMP(M,N) ,OUTP2(M,N) 

CHARACTER*40  TITLE 

COMMON/DAT/POWER , TEMP , MODEL , I START , I FI N , JSTART , JFIN 
COMMON  TO,TN,QM 
WRITE  ( I FN , 300 )  TITLE 
30(J  FORMAT  (1X.40A) 

DO  100  K=l, LINES 
I l=COLM* ( K-l ) +1 
I2-I1+COLM-1 
IF  (K.EQ. LINES)  I2=IFIN 
WRITE  (1FN.400)  (1,1=11,12) 

400  FORMAL  IX,  12  110) 

DO  200  J= JSTART, JFIN 

WRITE  (IFN.500)  J , (OUTP2( I , J ) , 1=11 , 12 ) 

500  FORMAT (IX, 13, IX, 12  E10.3) 

200  CONTINUE 

100  CONTINUE 

RETURN 
END 


APPENDIX  D 


THERMAL  ANALYSIS  FOR  ATHEROMATOUS  PLAQUE  ABLATION 

The  priaary  Beans  of  piaque  ablation  using  the  aicrowave 
angioplasty  technique  is  through  a  phase  change  of  the  plaque.  As  the 
fatty  tissue  is  heated  to  teaperatures  above  100 *C,  the  plaque  starts 
to  aelt  or  vaporize.  During  a  phase  change  the  teaperature  reaains 
constant  while  the  substance  is  in  equilibriua  at  its  Belting  or 
vaporization  point.  When  the  latent  heat  of  fusion,  vaporization  or 
subliaation,  is  reached,  soae  of  the  plaque  Bass  is  lost  due  to 
Belting  or  vaporization.  Assuaing  all  phase  changes  in  the  plaque 
takes  place  next  to  the  blood  and  plaque  boundary,  the  latent  heat  is 
absorbed  by  the  blood  (if  its  teaperature  is  less  than  that  of  the 
plaque) . 

The  bioheat  governing  equation  can  be  aodified  to  include  the 
effects  of  atheroaatous  plaque  ablation  by  eliainating  the  transient 
tera  and  replacing  it  with  a  latent  heat  accuaulation  tera. 

Eliaination  of  the  transient  tera  during  ablation  siaulates  the 
teaperature  being  constant  until  latent  heat  level  required  to  achieve 
plaque  ablation  is  reached. 

=  V(kVT)  +  qb  ♦  q,  (D-l) 

h..f?  =  *<«*>  ♦  Qb  ♦  0. 


(D-2) 


95 


The  differential  equation  can  be  rewritten  into  an  explicit 
finite  difference  equation  as  previously  discussed  in  Chapter  2. 
Again,  the  ablation  equations  are  derived  by  performing  an  energy 
balance  on  each  finite  difference  cell  and  using  the  physical 
formulation  method.  As  an  example,  the  explicit  finite  difference 
energy  balance  equation  for  the  general  convective  boundary  method  is 
shown  in  Figure  45.  Other  finite  difference  equations  at  the  various 
boundary  and  interior  nodes  are  similarly  derived. 

As  the  plaque  is  removed  from  the  model,  energy  is  conserved  by 
dumping  the  heat  into  the  blood.  This  increase  in  energy  is 
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Figure  45.  Derivation  of  the  plaque  ablation  explicit  finite 
difference  equation  for  the  convective  boundary  node  (Subroutine 
ABLCONV ) 
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represented  by  an  increase  in  blood  temperatures.  Some  heat  say  be 
conducted  further  into  the  plague  layer;  however,  this  amount  is 
assumed  to  be  relatively  small  since  the  melted  or  vaporized  plaque 
will  be  carried  away  be  the  blood  flow.  In  cases  where  only  one-half 
to  three-fourths  of  the  plaque  in  a  model  node  exceeds  the  critical 
latent  heat,  the  temperatures  will  be  allowed  to  increase  until  the 
entire  node  can  be  eliminated. 

The  type  of  phase  change  occurring  in  the  plaque  is  dependent 
upon  the  rate  of  heating.  The  critical  rate  of  heating  to  distinguish 
between  melting  and  heating  has  not  been  studied.  Estimates  for  the 
latent  heat  of  fusion  and  vaporization  were  made  using  the  thermal 
properties  for  fat  or  lipids.  As  future  experimental  studies  of 
atheromatous  plaques  are  completed,  provisions  have  been  made  in  the 
FDATA  computer  code  to  alter  the  default  values  for  fat. 

For  the  one  dimensional  fluid  flow  phaBe  changes  in  the  blood  can 
be  modelled  in  a  similar  model  as  shown  in  Figure  46.  This  subroutine 
has  not  yet  been  added  to  the  ablation  program  because  high  heat 
generation  within  the  blood  was  not  initially  anticipated  until  the 
results  from  the  FDTD  program  were  obtained. 

Like  the  ablation  subroutine,  blood  temperatures  are  held 
constant  until  the  heat  stored  in  the  cell  is  greater  or  equal  to  the 
latent  heat  of  vaporization.  The  blood  temperatures  are  then  allowed 
to  increase  as  more  heat  is  added  to  the  blood  by  microwave  heat 
generation  or  plaque  ablation. 
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AM  [■cp(T,_i-Ti)-«cp(Tim-Ti)] 
l.cc  =  h*Jt  =  — - ST1 - 1  +  2Trr0h.(T.-T,)  ♦ 

2*^(1,-!,)  +  Q.[  w<r0z-rj2)] 
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AM  =  pAV  =  P*(rp2-r.2)Ax 
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Figure  46.  Derivation  of  the  explicit  finite  difference  equation  for 
the  blood  vaporization. 


APPENDIX  E 


GLOSSARY  OF  MEDICAL  TERMS 
(Miller  and  Keane,  1987) 


ablation:  removal,  especially  by  cutting.  In  the  case  of  the 
microwave  angioplasty,  plaque  renoval  is  accomplished  by  heating 
tissues  to  the  point  where  a  phase  change  such  as  vaporization  or 
■elting  occurs. 

adventitia:  the  outer  coat  of  an  organ  or  structure,  especially  the 
outer  coat  of  an  artery. 

angina  pectoris:  acute  pain  in  the  chest  caused  by  interference  with 
the  supply  of  oxygen  to  the  heart. 

angioplasty:  plastic  repair  of  blood  vessels  or  lymphatic  channels. 

aorta:  the  largest  artery  arising  from  the  left  ventricle  of  the 
heart. 

arteriosclerosis:  thickening  and  loss  of  elasticity  of  the  coats  of 
the  arteries,  with  inflammatory  changes  caused  by  fibrous  and  mineral 
deposits  in  the  middle  layer  of  the  artery  wall. 

artery:  a  vessel  through  which  the  blood  passes  away  from  the  heart  to 
various  parts  of  the  body.  The  wall  of  an  artery  consist  typically  of 
an  outer  coat  (tunica  adventitial),  a  middle  coat  (tunica  media)  and 
inner  coat  (tunica  intima). 

atheroma:  an  abnormal  mass  of  fatty  or  lipid  material  existing  as  a 
discrete  deposit  in  an  arterial  wall. 

atheromatous  plaque:  a  deposit  of  predominantly  fatty  material  in  the 
lining  of  blood  vessels  occurring  in  atherosclerosis. 

atherosclerosis:  a  condition  characterized  by  degeneration  and 
hardening  of  the  walls  of  the  arteries  and  sometimes  the  valves  of  the 
heart,  related  especially  to  thickening  of  the  tunica  intima  (intimal 
layer)  due  to  an  accumulation  of  atheromatous  plaques. 

blood:  the  fluid  that  circulates  through  the  heart,  arteries  and 
veins.  It  is  the  chief  means  of  transport  within  the  body. 

catheter:  a  tubular  instrument  of  rubber,  plastic,  metal  or  other 
material,  normally  used  for  draining  or  injecting  fluids  through  a  body 
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passage.  For  coronary  angioplasties  they  serve  as  a  Beans  of 
delivering  the  plaque  removal  device  to  the  site  of  occlusion. 

cerebral  vascular  accident:  a  disorder  of  the  blood  vessels  of  the 
brain  resulting  in  impaired  blood  supply  to  parts  of  the  brain;  also 
called  a  stroke. 

cholesterol:  the  principle  aniaal  sterol,  occurring  in  faintly  yellow, 
pearly  leaflets  or  granules  in  all  aniaal  tissues. 

coronary:  a  tera  referring  to  the  vessels,  ligaaents,  etc. 

coronary  arteries:  two  large  arteries  that  branch  froa  the  ascending 
portion  of  the  aorta  and  supply  all  of  the  heart  auscle  with  blood. 

coronary  throabosis:  foraation  of  a  clot  in  a  coronary  artery. 

feaoral:  pertaining  to  the  feaur  or  to  the  thigh. 

heaoglobin:  the  oxygen-carrying  pigaent  of  the  blood,  the  principle 
protein  of  the  erythrocyte,  or  red  blood  cell. 

heaolysis:  rupture  of  erythrocytes  (red  blood  cell)  with  the  release 
of  heaoglobin  into  the  plasaa. 

heaatocrit:  the  volume  percentage  of  red  blood  cells  in  whole  blood. 

intiaa:  the  innermost  coat  of  a  blood  vessel;  also  called  tunica 
intima. 

in  vitro:  referring  to  studies  performed  under  artificial  conditions 
in  the  laboratory  on  tissues  removed  from  a  living  organism. 

in  vivo:  referring  to  studies  performed  on  tissues  and  organs  not 
removed  froa  the  living  organism. 

ischemia:  deficiency  of  blood  in  a  part,  due  to  functional 
constriction  or  actual  obstruction  of  a  blood  vessel. 

lesions:  any  daaage  to  a  tissue. 

luaen:  the  cavity  or  channel  within  a  tube  or  tubular  organ,  as  a 
blood  vessel  or  the  intestine. 

aedia:  aiddle,  especially  the  aiddle  coat  of  a  blood  vessel,  or 
tunical  aedia. 

aicrowave:  a  wave  typical  of  electromagnetic  radiation  between  far 
infrared  and  radiowaves. 

myocardial:  pertaining  to  the  auscular  tissue  of  the  heart. 

ayocardial  infarction:  foraation  of  an  infarct  (a  localized  area  of 
ischeaic  necrosis  produced  by  occlusion  of  the  arterial  supply  or  the 
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venous  drainage  of  the  part)  in  the  heart  Muscle,  due  to  interruption 
of  the  blood  supply  to  the  area. 

necrosis:  death  of  a  cell  or  group  of  cells  as  the  result  of  disease 
or  injury. 

occlusion:  the  act  of  closure  or  state  of  being  closed.  In  the  case 
of  coronary  occlusion,  this  tern  refers  to  the  obstruction  blood  flow 
through  an  artery  of  the  heart. 

percutaneous:  perforned  through  the  skin. 

perfusion:  the  act  of  pouring  through  or  into;  especially  the  passage 
of  a  fluid  through  the  vessels  of  a  specific  organ  or  body  part. 

recanalisation:  dilation  of  a  previously  occluded  blood  vessel  to 
restore  nornal  blood  flow. 

restenosis:  recurrent  stenosis,  especially  of  a  cardiac  valve  after 
surgical  correction  of  the  prinary  condition. 

stenosis:  narrowing  or  contraction  of  a  body  passage  or  opening. 

sterol:  a  solid  alcohol  of  aninal  or  vegetable  origin,  having 
properties  like  the  fats. 

transluninal:  through  the  lumen  or  channel  of  a  blood  vessel. 

tissue:  a  group  or  layer  of  similarly  specialized  cells  that  together 
perform  certain  special  functions. 

tunica:  a  covering  or  coat;  used  in  anatomic  nomenclature  to  designate 
a  membraneous  covering  of  an  organ  or  a  distinct  layer  of  the  wall  of  a 
hollow  structure,  as  a  blood  vessel. 

vasospasm:  a  sudden  involuntary  contraction  of  a  blood  vessel. 

vein:  a  vessel  through  which  blood  passes  from  various  organs  or  parts 
back  to  the  heart,  in  the  systemic  circulation  carrying  blood  that  has 
given  up  most  of  its  oxygen. 
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